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Abstract 

Magnetothermal instabilities are one of the peculiar phenomena of interest in conventional type-II, as well 
as in high-T c superconductors. In the present paper we attempt to analyze the nature and origin of the 
magnetothermal instabilities of the critical state and flux jumps phenomena in superconductors in the light 
of recent theoretical and experimental results. 
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Introduction 

The magnetic flux penetration into a type-II superconductor 
occurs in the form of quantized vortices. Each of these vor- 
tices has a normal core, which presents a long cylinder with a 
radius comparable to the superconducting coherence length £. 
Undamping superconducting current flows around this cylinder. 
The vortex current seizes an area of radius of order of penetra- 
tion depth A. Each vortex carries one magnetic flux quantum, 
which determined by the equality (f>o = 2 4- 10 7 Gs ■ sm 2 . The 
penetration of the vortex into the superconductor becomes en- 
ergetically favorable at the field H c \. The vortices are located 
one after another on the distance A forming correct triangular 
lattice in the superconductor at the interval H c i <C H <C H C 2; 
where H c \ and H C 2 are the lower and upper critical magnetic 
fields. With increasing magnetic field the period of lattices de- 
creases and the vortex density increases. At the field H C2 the 
vortex density becomes so great that the period of lattices be- 
comes an order of coherence length £. This means that magnetic 
flux completely penetrates into the superconductor and second 
type phase transition to the normal state occurs. In the interval 
of fields H c i <C H <C H C 2 superconductor is in the mixed state. 
In the presence of different types of defects or pinning centers 
in the superconductor sample, the vortices may be attached to 
such defects. A nature of interaction between the vortices and 
the structural defects is determined by the pinning force Fp. If 
transport current with the density j is passed through supercon- 
ductor, the interaction of the current with vortex lines leads to 
the emergence of the Lorentz force Fl, acting on each one of the 
vortices. Under the effect of the Lorentz force Fl the viscous 
flux flow of vortices begin to move. The viscous magnetic flux 
flow in accordance with electromagnetic induction creates a vor- 
tex electric field E. This means that energy dissipation occurs, 
an electric resistance appears and the superconductor undergoes 
a transition to the resistive or to the normal state. Propagating 
magnetic flux causes Joule heating, giving rise to global and/or 
micro flux avalanches in the critical state of type-II supercon- 
ductors. Thus flux jumps results in a large-scale flux avalanches 
in a superconductor and their origin are related to the magne- 
tothermal instabilities. 

Magnetothermal instabilities of the critical state and flux 
jump phenomenon in hard superconductors with high values of 
the critical current density and the critical magnetic fields have 
been investigated since 1960's in a classical works of Bean and 
et al. [1-3]. A detailed qualitative analysis of the magnetic in- 



stabilities in type-II superconductors was given by Wipf [4] . The 
magnetic flux jumps were studied in Refs. [5, 6] for a linear 
voltage-current characteristics of superconductor in the frame- 
work of adiabatic approximation i.e., assuming that the thermal 
diffusion, is much smaller than the magnetic flux diffusion. The 
criterion for the stability of the critical state in the case of dy- 
namic approximation was obtained by Kremlev [7] . The general 
concept of the magnetic instabilities in type-II superconductors 
was developed in literature [8] . The thermal and electromagnetic 
processes, whose development leads to the flux jump, have been 
investigated in detail and the stability criterion have been found 
in the framework of adiabatic and dynamic approximations in 
the viscous flux flow regime by Mints and Rakhmanov [9]. The 
detailed theoretical analyze of the flux jumping in the flux creep 
regime where the current-voltage characteristics of type-II su- 
perconductors is a nonlinear have been carried out recently by 
Mints [10] and by Mints and Brandt [11]. 

The dynamics of magnetothermal instabilities has been re- 
cently extensively studied in hard type-II as well as high T c super- 
conductors. And several contrasting models have been proposed 
for the detailed description for the nature of their origin and 
dynamics. The dynamics and the nature of observed in recent 
magnetization measurements of various techniques flux jumps 
can be explained in the scenario of magnetothermal instability 
theory [1]. According to this theory moving flux lines increase 
the local temperature T and therefore reduce the critical cur- 
rent density j c , which triggers further flux motion. This positive 
feedback may result in a magnetic flux avalanches and thermal 
instabilities in the superconductor sample. It is notice able that a 
recent theoretical and experimental investigations show that the 
dynamics of these flux jump instabilities strongly depends on 
the local temperature change, the rate of field change, the criti- 
cal current density and its temperature and field derivatives, the 
ratio between thermal conductivity and heat capacity, the sam- 
ple geometry, cooling conditions, pinning and micro-structural 
properties of a considered superconducting material. 

Many experimental and theoretical investigations on mag- 
netothermal instabilities in low and high-T c superconducting 
materials with different shapes and geometries, including thin 
Nb, MgE>2 films, organic superconductors and layered com- 
pounds, single and textured crystals, rings and bulk samples have 
been recently reported. However, no systematic studies of this 
problem have been reported and there is little understanding 
as to how different internal and external parameters affect flux 
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jumping behaviors in a various superconducting samples. 

Interestingly that most magnetization measurements show 
that the flux jumps occur at lower magnitudes of temperature 
and magnetic fields in a single superconductor material. Some 
experimental results show that flux jumps occur at temperatures 
and magnetic fields larger than the full penetration field H p , but 
lower than the H c \ . For instance, magnetothermal instabilities in 
organic type-II superconductors have been observed in the form 
of abrupt flux jumps at extremely low temperatures by Mola et 
al. [12] by means a torque magnetometer. Based on analysis of 
the temperature dependence of the jumps they found that the 
amplitude of the observed flux jumps increases with tempera- 
ture as A ~ T 3//2 , which consistent with accepted models. They 
have showed that as the temperature is increased the flux jumps 
are observed less frequently, and with greater magnitude. Thus, 
the authors believed that the sudden cessation of the flux jumps 
above a characteristic field Bj, which decreases upon increasing 
the temperature. It was also found that the average amplitude 
of the flux jumps is an increasing function of angle having an 
A(Q) « 1 — cosh(O) dependence. 

Similar temperature and field dependencies of the flux jump 
instabilities have been observed by Radovan and Zieve [13] in 
their local Hall probe measurements. They found both large and 
small flux jumps at different temperatures in a type-II thin Pb 
film superconductors. The flux jumps occur only at very low 
applied fields. A large jumps observed only at relatively high 
temperatures. It has been shown by these authors that at suffi- 
ciently low temperatures the jump size decreases roughly as T 3 , 
approaching a finite value at zero temperature. They found a 
power law distributions of avalanche size at temperatures T = 0.3 
K and T = 4.3 K with corresponding exponents of 2.01 and 1.09, 
respectively. Thus the authors argued that an interplay between 
vortex density and the microstructure is at the origin of the flux 
instabilities. 

The strong temperature dependence of the magnitude and 
frequency of flux jumps have been observed by Nowak and his 
co-workers [14] in a ring-shaped Nb thin films using a Hall probes. 
The magnetization of the film was measured at T = 1.4 4- 10 K 
as a function of temperature and applied field, perpendicular to 
the plane of the sample. At T ~ 3.1 K has been detected a 
crossover from a broad to narrow of the size distribution of in- 
ternal avalanches in the rings. The authors argued that the small 
value of specific heat at low temperatures is at the origin of these 
jumps in magnetization curves. So, these low-temperature huge 
avalanches in the rings may be qualitatively understood in a 
thermal instability scenario [9], according which at low temper- 
atures the vortex motion can result in a local temperature rise. 
The temperature rise leads to a decrease in the critical current 
density and therefore magnetothermal instability is produced. 
The stability parameter [8, 9] numerically has been calculated 
by using the measured temperature dependencies of the critical 
current density j c and heat capacity v ~ T . The authors found 
that at highest temperatures thermal effects are negligible and a 
largest in magnitude heat capacity serves as a stabilizing factor 
of the Bean's critical state. 

A such type of flux jumps in the magnetization curves at 
low temperatures for a single YBaCuO crystals with a different 
transition temperatures T c has been found by Khene and Bar- 
bara [15]. It was shown that the first flux jump field Bj depends 
not only on critical current density j c and specific heat v, but also 



on the field sweep rate and the transition temperature T c of the 
sample. The results of their magnetization measurements show 
that the field for the first flux jump Bj increases as temperature 
increases and the flux jumps disappear at higher temperatures. 
The results obtained in this work suggest a better thermal and 
magnetic stability for considered single crystals with low tran- 
sition temperature T c . The effect of magnetic field sweep rate 
on the behavior of flux jumps has also been investigated. The 
field Bj decreases with increasing magnetic field sweep rate B e . 
Interestingly that these results are very similar to those obtained 
in conventional II-type superconductors [16]. 

A more detailed the temperature and magnetic field sweep 
rate dependencies of flux jumping has been analyzed by Nabi- 
alek et al. [17] at temperature of 4.2 K in poly crystalline single 
crystals by means of magnetization, screening and torque mea- 
surements. It has been shown that as the sweep rate decreases 
the value of the first flux jump field Bj increases rapidly and 
tends to saturate at higher sweep rates. Similar strong sweep 
rate dependence on the first flux jump field Bj early has been 
observed in other experiments of McHenry [18] and Guillot [19]. 
Nabialek et al. [17], also investigated the temperature depen- 
dence of flux jumping, as well as the influence of flux creep and 
demagnetizing effects on the critical state stability. 

The behavior of flux jumps with different amplitudes in su- 
perconductors have studied Vanacken et al. [20] using pulsed 
field magnetization measurements. The time dependence of the 
flux jumps was presented. A periodical jumps were observed for 
both the positive and negative field polarities. Surprisingly that 
the sharp flux jumps occur at lowering magnitude of the field 
for all temperatures below T = 30 K. These jumps are expected 
to be equidistant Bj £s 0.3 T. As has been shown that the flux 
jumps strongly depend on the magnetic field sweep rate and the 
magneto-thermal history of the sample, on the temperature and 
the critical current density. Such dependence of occurrence of 
flux jump field on the sweep rate of external magnetic field and 
the temperature for the melt-textured superconductor materials 
has been studied by Fuchs et al. [21] at temperatures T > 20 
K. The observed magnetic instabilities were discussed within the 
framework of existing theoretical results. 

Recently Zhao et al. [22] have observed many multiple small 
and irregular flux jump instabilities in MgB2 thin films at tem- 
perature below 2 K and applied field of 1.3 T by means a SQUID 
magnetometer measurements. The magnitude of these instabil- 
ities are much larger in low fields than in high fields. With in- 
creasing of temperature the small instabilities evolved into some 
larger ones and disappear completely at temperature T = 14 K. 
The authors believed that there are many places for the avalanche 
to occur due to a high density of small defects formed during the 
preparation process of the thin films leading to much stronger 
critical current densities. As pointed authors thermal diffusion 
is much easier in thin films due to their very small thickness 
and large surface area exposed to the environment and conse- 
quently, in thin films each avalanche is small in magnitude but 
the number of avalanches can be huge. The authors argued that 
this simple picture may give an explanation to many small flux 
jumps observed in superconducting thin films. They reported 
the observation of the suppression of the critical current density 
at low temperatures due to many small jumps in thin films. 

Muller and Andrikidis [23] using a SQUID magnetometer has 
measured the magnetization loop for a melt-textured high T c su- 
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perconductors at different temperature and magnetic field inter- 
vals. The first flux jump occur at temperature T=3.0 K in the 
magnetization loop, when the applied field is perpendicular to 
the crystallographic c axis of the sample. At temperature T=3.5 
K solitary jumps were found and above 3.9 K no flux jumps oc- 
cur in the magnetization curve. With increasing temperature 
the first jump field increases and at higher temperatures (above 
T=3.05 K) flux jump suddenly disappears. The measured val- 
ues of Bj has been compared with existing theoretical results [8, 
9]. In this work the stability criterion for flux jumps has been 
analytically calculated using Kim- Anderson critical state model. 

Chabanenko at al. [24] have studied magnetothermal insta- 
bilities and giant flux jumps, both theoretically and experimen- 
tally in the framework of adiabatic approximation using vari- 
ous dependencies of the critical current density on the magnetic 
field. In particular, they have numerically calculated magnetiza- 
tion and magnetostriction loops with flux jumps employing the 
Kim- Anderson critical state model and exponential model for the 
dependence j c (H). On the basis of these theoretical calculations 
it has been constructed a H-T phase diagram of flux jump insta- 
bilities. The authors showed that the maximum value of tem- 
perature after the flux jump strongly depends on the magnetic 
field sweep rate and micro-structural properties of the sample. 

Vasilev et al. [25] experimentally studied the time of duration 
of the flux jumps and amount of the magnetic flux entering into a 
disc type of superconductor sample as a function of temperature 
and external magnetic field. In another paper Chabanenko et al. 
[26] have studied the structure of the flux jumps in different su- 
perconducting samples, such as Nb-Ti samples, poly crystalline 
Nb plates and melt-textured YBaCuO slabs on the basis of a 
temporal dependence of the surface magnetic induction using a 
miniature Hall probe sensors. The authors have observed very 
interesting oscillating phenomena originated by the thermomag- 
netic avalanches in the vortex matter for both low and high-T c 
superconductor samples. These oscillation can be qualitatively 
interpreted in terms of the theoretical model which takes into 
account the existence of a definite value of the effective vortex 
mass, i.e the inertial properties of the vortex system [27]. In this 
work a details of three stages of the thermomagnetic instability 
development was also proposed. 

The influence of the time-dependent external conditions on 
the stability threshold of the critical state for hard type-II 
superconductors has been studied theoretically by Mints and 
Rakhmanov [28] . They found a stability criterion of the critical 
state as a function of the rate of variation of the external mag- 
netic field and the external temperature in hard superconductors. 
The authors argued that the existence of a background electric 
field induced by a variable external magnetic field is essential for 
observing a temperature and electric field oscillations in a super- 
conductor sample. Similar magnetothermal oscillations in the 
form of flux jumps have been observed by Kumm et al. [29] in a 
melt-textured crystals in external varying magnetic field. The in- 
terval AB e between consecutive jumps increases with sweep rate 
B e and decreases with temperature. Every jump is preceded 
by oscillations of the induced voltage and the sample tempera- 
ture. A similar oscillations preceding the flux jump instabilities 
have been studied by numerous researchers in earlier both ex- 
perimentally by Shimamoto [30] and theoretically by Maksimov 
and Mints [31]. Earlier Zebouni et al. [32] have detected that 
the period of oscillations increases strongly with increasing tem- 



perature. The oscillations disappear sharply if the sweep was 
stopped. This effect due to the simultaneous and sharp drop of 
the sample temperature to the bath temperature. Such a series 
of heating pulses of very large amplitude appeared at discrete 
reproducible values of B. The pulses were observed in decreas- 
ing field, but at relatively lower values of B. Simultaneously 
with each of these heat pulses, sudden penetration of magnetic 
flux was observed. Therefore the pulses related to discontinuous 
steps in the magnetic flux penetration - flux jumps. A tenta- 
tive identification of these fluctuations with one of the collective 
modes of vibration of the vortices predicted by de Gennes and 
Martison [33]. Such limited flux jump instabilities, initiated by a 
varying external magnetic field perturbation has been observed 
in ceramic high-T c superconductors by Bodi et al. [34]. An 
extensive study of the magnetothermal oscillations and limited 
flux jumps in a high-T c superconductor samples was given by 
Legrand et al. [35]. Later Legrand et al. [36] have proposed 
that the magneto-thermal oscillations in the Bean's critical state 
may arise as coupled oscillations of small perturbations of the 
temperature and the electric field. They both theoretically and 
experimentally studied the dependence of the frequency of the 
oscillations on the magnetic sweep rate and the temperature for 
a granular superconductors. The experimental results for mag- 
netothermal oscillations in the temperature T = 3 ~ 7 K and 
magnetic field sweep rate B e — 10 -f - 45 G/s ranges show that at 
low values of the applied field B e (t) there is slow temperature in- 
crease due to the small vortex avalanches establishing the critical 
state. Above certain threshold magnetic field, will appear tem- 
perature oscillations with the period t in the range t ~ 10 -t- 70 
s and its amplitude increasing in time. A flux jump occurs close 
to the Bean field accompanied by a temperature rise up to about 
12 K with a characteristic time of the order of 1 s. The frequency 
of these oscillations is proportional to the magnetic field sweep 
rate B e . 

Recently, it has been shown [37] that the local Joule heat- 
ing due to planar defects or low-angle grain boundaries can 
cause thermal instabilities in a high-Tc coated conductors. It 
was shown that the thermal instabilities are more pronounced 
at lower temperatures and magnetic fields or under poorer gas 
cooling. The author have obtained the stability criterion for 
a planar defect in a thin film for a power-law voltage-current 
characteristics of the sample on the basis of calculation a steady- 
state heat balance equation. It has been pointed by Gurevich 
[37], also that the highly nonlinear voltage-current characteris- 
tics may cause strong disturbances of the electric field and dis- 
sipation with spatial size much larger than the defect size in the 
sample. Thus the resulting local hot spots near defects can trig- 
ger thermal instability in thin films. At sufficiently low electric 
fields E the voltage-current characteristics of type-II supercon- 
ductors is highly nonlinear and differential conductivity a is a 
function of the electric field E [38, 39]. The dependence a of E 
at small values of the electric field E may considerably affect the 
critical state stability and the occurrence of the flux jumps in the 
sample [39]. As recently has been shown by Mints [10] that the 
nonlinear differential conductivity a significantly affects the flux 
jumping process and therefore the first flux jump field Bj in the 
flux creep regime. In particular, the author derived a criterion 
for thermomagnetic flux-jump instability for the case of trans- 
verse geometry of thin films taking into account the nonlinear 
dependence of the background electric field on sample conduc- 
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tivity in the framework of Bean's critical state model. It was 
shown that the magnetic field sweep rate affects the instability. 
This effect is connected with the non-linear electric field depen- 
dence of j(E) in the flux creep regime. The first flux jump field 
Bj decreases with increasing sweep rate field B e . As flux creep 
is relatively strong in high-T c superconductors, one may expect 
it to have a significant influence on the critical state stability in 
these materials. 

A detailed measurements have been performed by means a 
torque and SQUID magnetometer during magnetic field sweep 
rate by Monier and Fruchter [40] in a single crystal of the or- 
ganic superconductor. They found an expression for the first flux 
jump field Bj for a cylinder sample assuming that the temper- 
ature in the sample is uniform. The obtained criterion for flux 
jumps for a nonlinear voltage-current characteristic in the dy- 
namic approximation in a good agreement with a recent results 
of Mints [10]. It was shown that magnet othermal instabilities 
may be suppressed by an increase of heat transfer coefficient or 
the use of thin samples. The activation energy, the screening 
current density and its temperature derivative has been found 
from the dynamic relaxation measurements. 

In the present paper, we shall investigate the magnetother- 
mal flux jump instabilities of the critical state in conventional 
and high-r c superconductors. We determine the flux jump field 
Bj and critical state stability criterion within the framework of 
Bean's model. To determine the flux jump field an analytical 
simulations of coupled equations for the magnetic B(r, t) and 
electric E(r, t) field inductions and temperature T(r, t) has been 
performed. The field of the first flux jumps Bj is calculated ana- 
lytically by using the dynamic and the adiabatic approximations. 
The qualitative analysis of the magnetic flux jumps instabilities 
for superconductors will be provided. 

§1. Flux jump instabilities 

Let us suppose that at the initial moment of time, the mag- 
netic field in the sample is uniform and equal to Bo , after which 
the external field rises to some value B e . With an increase of 
the external field B e , the vortices penetrates the sample. The 
viscous flow of vortices inside the sample causes an electric field 
E which generates persistent currents near the sample surface, 
having a density j ~ j c and power dissipation. This dissipation 
heats the sample which leads to further changes in the magnetic 
flux density and further increase the sample temperature locally 
decrease in the critical current density j c . With a decrease j c , 
the vortex lines penetrates the sample more deeply, more heat 
is released and new dissipation occurs thus resulting in a new 
increase of temperature T. Under certain conditions such a dis- 
sipative flux motion can cause magnetothermal instabilities. 

Mathematical problem of theoretical study the dynamics of 
thermal and electromagnetic perturbations in a superconductor 
sample in the flux creep regime can be formulated on the basis of 
a system nonlinear diffusion-like equations for the thermal and 
electromagnetic field perturbations with account nonlinear rela- 
tionship between the field and current in superconductor sample. 
The distribution of the magnetic flux density B and the trans- 
port current density j inside a superconductor is described by 
the equation 

rotB = — j. (1) 



When the penetrated magnetic flux changes with time, an elec- 
tric field E is generated inside the sample according to Faraday's 
law 

,2, 

The temperature distribution in superconductor is governed by 
the heat conduction diffusion equation 



v(T) a ±=V[n(T)VT]+jE, (3) 



Here v = v(T) and k = k(T) are the specific heat and ther- 
mal conductivity, respectively. The above equations should be 
supplemented by a current-voltage characteristics of supercon- 
ductors, which has the form 



j=j c (T,B)+j(E). (4) 



In general, the critical current density depends, on both the lo- 
cal magnetic field B and local temperature T, thus, the magnetic 
flux profile is determined by equation j c — j c (B,T). CP. Bean 
[2] has proposed the critical state model which successfully used 
to describe magnetic and transport properties of hard supercon- 
ductors. According to this model the current density j is equal 
to the critical current density j c and independent of magnetic 
field. In order to obtain analytical results of a set Eqs. (l)-(4), 
we suggest that j c is independent on magnetic field B and use 
the Bean critical state model j c = j c (B e ,T). We suppose that 
the critical current density j c is linearly dependent on the local 
temperature 



jc(T) = jo — a(T — To). (5) 

where a = jo/(T c — To); jo is the equilibrium current density; 
To and T c are the equilibrium and critical temperatures of the 
sample, respectively. The characteristic field dependence of j(E) 
in the region of sufficiently strong electric field (E > Ef) can 
be approximated by the piece-wise linear function j = a/E; 
<jf — r)c 2 1 B$o = u n H C2 /B is the effective conductivity in the 
flux flow regime; r) is the viscous coefficient, $o = irhc/2e is 
the magnetic flux quantum, o n is the differential conductivity 
in the normal state, H C 2 is the upper critical magnetic field, Ef 
is the boundary of the linear area in the voltage-current charac- 
teristics of the sample [41]. Notice, that in the flux flow regime 
the differential conductivity 07 is approximately constant, i.e., 
independent on the electric field E. 
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The boundary conditions for the magnetic field induction are 




Fig.l The geometry of the problem 

Let us formulate a differential equations governing the dynamics 
of small temperature and electromagnetic field perturbation in a 
superconductor sample. We study the evolution of thermal and 
electromagnetic penetration process in a simple geometry - su- 
perconducting semi-infinitive sample x > (Fig 1.). We assume 
that the external magnetic field induction B e is parallel to the z- 
axis and the magnetic field sweep rate B e is constant. When the 
magnetic field with the flux density B e is applied in the direction 
of the z-axis, the transport current j(x,i) and the electric field 
E(x, t) are induced inside the slab along the y-axis. For this ge- 
ometry the spatial and temporal evolution of small thermal T(x, 
t) and electromagnetic field E(x, t) perturbations are described 
by the thermal diffusion equation coupled to Maxwell's equations 



d 2 E 
dx 2 



'lit 



, d 2 T . ^ 
k dx^ + ^ 



(6) 
(7) 



dj dE dj c dT 
_dE~dt ~ dTltt 
It should be noted that the nonlinear diffusion-type equations 
(6) and (7), totally determine the problem of the space-time dis- 
tribution of the temperature and electromagnetic field profiles in 
the flux flow regime with a linear current-voltage characteristics 
in the semi-infinite sample. 

The system of differential equations (6) and (7) may be inte- 
grated analytically subject to appropriate initial and boundary 
conditions for the thermal and electromagnetic small perturba- 
tions. We present the thermal boundary conditions as 

dT{0,t) 
dx 

where w is the heat transport coefficient. Let us assume that the 
magnetic field perturbation is equal to zero at the sample surface 
and according to relation (2), we obtain the first electrodynamic 
boundary condition 

dE(0,t) 



-w(T-To), T(L,t) = T 



(8) 



dx 



0. 



(9) 



The second boundary condition for the electric field E(x, t) at 
the flux front x — L can be presented as 



E(L,t)=0, 



(10) 



dB(0,t) = B e , B(L,t) = 0, 



cB e 



where L = is the London penetration depth. 
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§1.1. Dynamic instability 

In this section we study the dynamics of flux-jump instability 
of during flux penetration into superconductor sample within the 
dynamical approximation [8]. As we know [8, 9] that a nature 
of the flux jumps depends on the competition between diffusive 
and dissipative processes through the dimensionless parameter 



471(7/ K 



Dt_ 
D m 



where Dt — n/v is the thermal diffusivity and D m = c 2 /4naf 
the magnetic diffusivity coefficients, respectively. Therefore the 
flux instability criterion is determined mainly by the relation of 
the magnetic D m and thermal Dt diffusion coefficients. Let us 
consider the case for the flux jumps corresponding to the limit- 
ing case t 1. Consequently, it can be assumed that the initial 
rapid heating stage of a flux jump takes place on the background 
of a "frozen-in" magnetic flux. Therefore, under this dynamic 
approximation, we obtain from (7) the relation between electric 
field E(x, t) and temperature T(x, t) perturbations in the follow- 
ing form 



dj 



: 



(11) 



dE T c - T ' 
We notice that the last relation between E(x, t) and temperature 
T(x, t) has been derived in the assumption that the decrease 
of the critical current density j c resulting from a temperature 
perturbation T(x, t) compensates with increase of the resistive 
current density j resulting from an electric field perturbation 
E(x,t), so a total current density remains constant [8]. Upon 
substituting the expression (11) into the equation (6) and ex- 
cluding the variable E(x, t) one can get the differential equation 
for the distribution of thermal perturbation, which can be con- 
veniently presented in the following dimensionless form 



d 2 o | de | SQ 

dz 2 dr 



0. (12) 
Here we introduced the following dimensionless variables 



e 



T 

and parameters 



t 

tk ' 



E e = 



ah 2 



tk 



8 = 



vL 

K 



o f E e 



x 
V 



The equation (12) can be easily solved by using a method of sep- 
aration of variables. Taking into account the thermal boundary 
conditions 

<Z0(O,t) 



6(l,r) 



dz 



0. 



(13) 



an explicit solution to equation (12) can be written as the fol- 
lowing form 



T(x,t) = (T c -T )exp 



4L 2 J 



(t - to) 



cos— x + T (14) 
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where to is a constant parameter. The distribution of the elec- 
tric field perturbation E(x, t) can be found by using relations 
(11) and (14) 



E{x, t) = E c exp 



IT 

4E 2 



(t — to) ) cosh —x (15) 



Using the obtained solutions (14) and (15) one can easily find 
the critical field Bj starting from which flux jump takes place 



Bi = 



16 



KOfJc 

a c 2 



-,1/2 



and the critical thickness d c of the sample 



aj c ' 



(16) 



(17) 



We can easily estimate Bj and d c using by a typical numerical 
values of a well known parameters of the sample which has been 
made in many publications (see, for example [8]). 

The dynamic approximation which we used here is justified 
if the thermal flux diffusion is much greater than the magnetic 



diffusion D t D m . The thermal D t 



1(T 



10 • cm / s and 



magnetic D m — 1 -j- 10 ■ cm 2 / s diffusion coefficients we esti- 
mated using by a numerical values of parameters. On the other 
hand, using the obtained numerical values of D t and D m it is 
easily estimate the thermal and magnetic diffusion times, taking 
into account that, both t K and t m are inversely proportional to 
the corresponding thermal and magnetic diffusion coefficients, 
respectively. Then one can be seen, within the time interval 
t K <C At <JC t m the results obtained provide for a highly accu- 
rate description of the dynamical evolution of temperature and 
electromagnetic perturbations in the superconductor sample. 

§1.2. Adiabatic instability 

In the adiabatic approximation we assume that rapid prop- 
agation of the flux line is accompanied by an adiabatic heating 
of the superconductor [8]. In this case the thermal conditions 
have a little effect on the occurrence of the flux jumps. There- 
fore thermal conductivity is negligible and corresponding term 
can be neglected in the heat diffusion equation. In this limiting 
case t « 1 the distribution of the electric field is described by 
the following nonlinear equation 



1 dr 2 



(21) 



where k is the constant parameter to be determined. The equa- 
tion (20) has an explicit solution 
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Let us rewrite the equation (21) in the form of 



which has an exact solution in the form 



fir) 



k y/2 



(22) 



(23) 



(24) 



where r p is the integrating constant and b is the free parameter. 
Combining above solutions (22) and (24) we get the following 
expression for the electric field and temperature distribution 



E(x,t) 



E e (L 2 — x 2 ) 
2 ^Ht-t^ 

V2 



T(x,t) = ^(L*-x>) 



ch 



< b(t-t p ) 
V2 



(25) 
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(26) 

According to (25) perturbation E(x, t) decreases exponentially 
with time At — y2/b. Motion of flux lines caused by E(x,t) 
perturbations leads to increase temperature in the region x < L 
of the sample. From (26) at t — t p ^> At we obtain the following 
expression for the maximum heating during the flux jump 
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(27) 



Numerical estimation gives AO m ~ 2 -j- 3 for the typical val- 
ues of parameters of the sample. Such amount of heating will 
be dissipated during the adiabatic flux jump instability which is 
determined by a well known stability parameter 
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The solution of (18) together with the boundary conditions (8) 
can be easily obtained analogously, as in a previous section, so 
we have 

e(z,r) = X(r)ip(z). 

after substituting the last relation into (18) one obtains the fol- 
lowing expressions for the functions <j>{x) and A(r) 



dx 2 



A more detailed derivation of the adiabatic stability criterion for 
Off flux jumps within the framework of Bean critical state model 

a has been given in [8]. To compare the theoretical value of j3 with 

experimental data, one must know the temperature dependence 
the heat capacity coefficient of the sample, which at low tem- 
peratures can be expressed as v{T) ~ vq(T /T ) n , where v is 
a numerical factor, n ~ 3.1. Using the expression for v(T) one 

(19) can estimate the value of j3 numerically for an existing of various 
functional dependencies of j c (T). It should be noted that due 
to high values of critical temperature and thus heat capacity in 
high T c superconductors, the experimental value of the first flux 
jump field /3 can be larger than in conventional type II supercon- 

(20) ductors. 
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§1.3. Flux creep 



the dependence j(E) has the form 



In general, the critical state stability threshold is determined 
by many external and internal factors, as by the sweep rate of the 
external magnetic field, the type of voltage-current characteris- 
tics, the critical current density and its magnetic field and tem- 
perature derivatives, the profile temperature and surface cooling 
conditions, the sample geometry and pinning properties of the 
considered sample. The obtained above results are valid in the 
flux flow regime, where voltage current-current characteristics of 
hard superconductor is described by linear dependence of j(E) 
at sufficiently large values of electric field [41]. As can be seen 
from the obtained results for the flux flow state the value of Bj 
depends mainly on the critical current density and the specific 
heat of the sample. The nonlinear part of the curve j(E) in the 
region of weak electric fields is associated with flux creep. In the 
flux creep regime the magnitude of flux jump field Bj strongly 
depends on variation of external parameters, in particular, on 
the magnetic field sweep rate [10, 11]. It can be shown [10, 11] 
that flux creep may sufficiently affect on the occurrence of the 
flux jump instabilities. In this section we shall discuss the effect 
of flux creep and the nonlinear current-voltage characteristics on 
the threshold of flux jumps, qualitatively. The current-voltage 
characteristics of type- II conventional as well as high-T c super- 
conductors in the flux creep regime is a highly nonlinear due to 
thermally activated magnetic flux motion. Thermally activated 
flux motion or flux creep problem in superconductor samples 
with various geometries and conditions has been recently exten- 
sively studied by many researchers [42-64]. According to Kim- 
Anderson theory [55, 56] the thermally activated flux motion is 
described by a well known expression 



v = v exp[-U/kT\, 



(29) 



Here Vo is the velocity of the thermally activated flux motion 
at zero temperature, U is the activation energy due to vortex 
pinning. The activation energy U = U(j, B, T) depends on tem- 
perature T, magnetic field induction B and current density j. 
The dependence of U on the current density j is extensively dis- 
cussed in the literature [57] . In particular, Burlachkov et al. [58] 
have analyzed the flux creep at different dependencies of the acti- 
vation energy U on field B and current density j. For the vortex 
glass and collective creep models the potential barriers highly 
nonlinear function of j [59]. A voltage-current characteristic of 
a type-II superconductor in the flux creep state is characterized 
by the power or exponential law increase of E with increasing j. 
For the linear current dependence of the potential barrier U (j) 
[55, 56], the dependence j(E) has the form 



: j c + ji In 
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eZ 



(30) 



where parameter ji determines the slope of the j-E curve and it 
is assumed ji <C j c . In this case the differential conductivity a 
is determined by the following expression (see, Refs. [10, 11]) 



dE 



E' 



(31) 



For the logarithmic current dependence of the potential barrier 
U(j) proposed by Zeldov et.all. [59] 



U(j) = U \n 
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3 J 



(32) 
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when the flux creep is determined by numerous spatial defects 
of the sample. J7o=const and E c is the crossover electric field. 
Here the parameter n = Uo/kT is a function of temperature T, 
magnetic field H and depends on the pinning regimes and can 
vary widely for various types of superconductors. In the case 
n = 1 the power-law relation (33) reduces to Ohm's law, de- 
scribing the normal or flux-flow regime. For infinitely large n, 
the equation describes the Bean critical state model j = j c [1]. 
When 1 < n < oo, this equation describes nonlinear flux creep 
[59] . In this case the differential conductivity a is determined by 
the following expression 



(34) 



a = d 3 = 3c 
dE nE' 

It is assumed, for simplicity, that the value of n temperature and 
magnetic-field independent. It should be noted that the nonlin- 
ear diffusion- type equations (6) and (7), completed by the flux 
creep equation (34), totally determine the problem of the space- 
time distribution of the temperature and electromagnetic field 
profiles in the flux creep regime with a nonlinear current-voltage 
characteristics in a semi-infinite superconductor sample. 

It should be noted that the investigation of the stability con- 
ditions in the flux creep regime is very difficult due to absence 
of solution of system equations (6)- (7) together with nonlinear 
j(E) dependence. However, in some limiting cases, it can be 
solved the problem if we take into account that the heating due 
to viscous flux motion is negligibly small and get some approxi- 
mate solution, describing the evolution of thermal and magnetic 
field diffusion in the creep regime. According to relation (34) 
the conductivity decreases with the increasing of electric field E, 
while a ~ E~ x and it strongly depends on the external magnetic 
field sweep rate E ~ B e . Therefore the stability criterion also 
strongly depends on the differential conductivity a [10]. 

§1.4. Dynamical approximation 

For the small thermal 5T(x, t) and electromagnetic field 
5E(x,t) 

5T = e(x)exp[yt], (35) 

SE = e(x)exp[jt]. (36) 

perturbations the system of differential equations (6) and (7) can 
be written in the following form 



vjO = k— — +j c e, 



dx 2 



d"e 
dx 2 



47T 



-7 



e 



(37) 
(38) 



3c Jc 

nE b e T c -T 

where 7 is the eigenvalue of the problem to be determined. It is 
clear that the rate 7 characterizes the time development of the 
instability. In the case, when Re7 > 0, small thermal and elec- 
tromagnetic perturbations increase and the stability margin cor- 
responds to the case when 7=0. It should be noted that the non- 
linear diffusion-type equations (37) and (38), totally determine 
the problem of the space-time distribution of the temperature 
and electromagnetic field profiles in the flux creep regime with a 
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nonlinear current- voltage characteristics in the semi-infinite sam- 
ple. 

As we have mentioned above, the differential conductivity 
cr(E), which determines the dynamics of the instability is high 
in the flux creep regime and the parameter r is high enough, 
also [10, 11]. It is clear that this picture for the flux jumps cor- 
responds to the limiting case r 3> 1. Consequently, it can be 
assumed that the initial rapid heating stage of a flux jump takes 
place on the background of a "frozen-in" magnetic flux. There- 
fore, under this dynamic approximation, we obtain from (38) 
the relation between electric field e(x, t) and temperature <d(x,t) 
perturbations in the following form 



Jc 



3c 



-9 = 



nE~ T c -T ~ " (39) 
Upon substituting the expression (39) into the equation (37) one 
can easily get a differential equation for the temperature distri- 
bution, which is conveniently presented in the following dimen- 
sionless form 



d 2 Q 



dp 2 



pQ = 0. (40) 
Here we introduced the following dimensionless variables 

1 1/3 



7 ■ 



x 



ab 
n E L 

Hi 



El ~ B e 



Thus, the condition of existence of a non-trivial solutions of equa- 
tion (40) allows to determine the spectrum of eigenvalues of 7 
and the instability threshold, accordingly. The equation (40) has 
an exact solution in terms of Airy functions given as the following 
form 

9(p) = ci(s)Ai(p) + c 2 (s)Bi(p). (41) 

where Ai(p) and Bi(p) are the Airy functions. Here constants of 
integration ci and C2 are determined from the thermal boundary 
conditions. Substituting the last solution (41) into the thermal 
boundary conditions we find that C2 = and O(p) = ci(s)Ai(p). 
Applying the second boundary condition we get an equation to 
determine the eigenvalues of the problem 



>2/3{a-n) = J-2/3( a n)- 

where a n are the zeros of the Bessel function, which grows with 
increasing n. For example, for n=l the stability criterion is pre- 
sented as 

01 = r 2/3 7 . (42) 

Using the value for the magnetic field penetration depth, we can 
easily obtain from (42) an expression for the threshold magnetic 
field Bj at which the flux jump occurs 



otiEl 



(43) 



Let us now estimate the threshold field for a typical values of 
parameters j c ^ W 9 A/m 2 , T c ~ T ~ 10 K, k ~ 10~ l W/Km, 

~ B e L, induced by the 
3 T/s is of the order of 



n=10. The background electric field El 



magnetic-field variation B e ~ 10 



10 



E L 



10" 



10~ J V/m for the value of L = 0.01 m. We can 



easily estimate that the threshold field has the value Bj 



l-=-3T. 



Experimentally, the background electric field El is created 
by the sweeping rate of the applied magnetic field B e . As can 



be seen from the relation (43) the threshold field Bj is decreased 
with the increasing of background electric field El- It is no- 
ticeable that the dependence of the flux-jump field Bj on the 
sweeping rate B e of the applied magnetic field have been veri- 
fied by a numerous experiments [15-21]. An intensive numerical 
analysis on the sweep rate dependence of the threshold field has 
been performed recently in [22] . Recent magnetization measure- 
ments [17] have shown that the value of the threshold field Bj 
decreases as the sweep rate B e increases. A theoretical investiga- 
tions on the dependence of threshold field on the varying external 
magnetic field has been performed by many researchers (see, for 
example Ref. [8]). Within the framework of the flux jump in- 
stability theory, a rapid variation of the applied magnetic field 
acts as the instability-driving perturbation, and that threshold 
field Bj should decrease with increasing the sweeping rate B e 
[8]. The numerical studies [22] have demonstrated that the flux 
jumps takes place when the sweep rate B e increases up to a cer- 
tain value, where the number of jumps increases with the sweep 
rate B e . As the sweep rate further increases, these simulation 
results show that the flux-jump field decreases and approaches a 
saturation value, which is fairly close to the experimental value 
of about 10 2 Oe [17]. However, as has been mentioned in [17], an 
experimental investigations on the dependence of the threshold 
field for the flux jump field Bj on the external magnetic field 
sweep rate B e is very little. Experimentally, can be observed a 
complex behavior of dependence of the threshold field Bj on the 
sweep rate B e . The results of experiments of Ref. [17] demon- 
strated that Bj is independent of the sweep rate in a defined 
range of temperatures. Thus, the near independence of Bj on 
the sweeping rate remains to be explained. In some conventional 
superconductors both the independence of Bj on the sweeping 
rate and its growth at a high sweeping rate [17-19] were detected. 
It has been suggested [9] that a nonuniform heating may be re- 
sponsible for such an effect. However, theoretical understanding 
of the thermomagnetic instabilities at such conditions is still lack- 
ing. Note, however, that some details of the local field behavior 
depend indeed on the sweeping rate, as, for example, the number 
and amplitude of the jumps. In [22], it has been demonstrated 
that at low values of the sweep rates B e the number of flux jumps 
decreases as sweep rate increases. At still high sweep rates the 
amplitude of flux jumps becomes independent of the sweep rate 
and saturates to the limit with further increasing sweep rate. 
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Fig. 2 and 3. The sweep rate dependence of the threshold field 

In Fig. 2 and 3 we have demonstrated the dependence of the 
threshold field Bj on the external magnetic field sweep rate at 

respectively. As can be seen, the 



B e 



-1/2 



and ~ B e 



-1/3 



value of Bj decreases as the sweep rate increases. As the sweep 
rate increases the value of Bj decreases and it tends to saturate 
at high sweep rates. Magnetic field dependence of the critical 
current density only slows down the decrease of the field Bj with 
increasing external magnetic field sweep rate [10, 11]. We note 
that for the case Kim- Anderson model [55], the absolute value 
of the exponent in the power formula decreases from 1/2 to 1/3, 

-1/3 



so Bi 



B e 



10 



Thus, the stability criterion for the flux jumps demonstrates 
extremely high sensitivity of the threshold field Bj on the values 
of the critical current density j c , thermal conductivity k, and 
it external magnetic field sweep rate B e . It follows from above 
criterion that the value of the threshold field Bj is inversely pro- 
portional to the square root of the magnetic-field sweeping rate 
Be . Therefore, with the increase of sweeping rate B e the thresh- 
old field Bj decreases, as can be seen from Figures 2 and 3. 



Since we have neglected the redistribution of heat in deriving 
(44), only, electrodynamic boundary conditions should be im- 
posed on this equation 



<1, 



de(0,r) 
dz 



(45) 



An explicit solution of (44) together with the boundary condi- 
tions (45) can be obtained by using the method of separation of 
variables. Looking for the solution of equation (44) in the form 



e( Z ,T) = \(r)iP{z). (46) 
we get the following expressions for a new variables 



dX 
rf7 



= -k\ L 
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dz 2 



(47) 



(48) 



By integrating equation (47) we easily obtain 



A(t) = (t, - r) 



-l/ 7 



(49) 



where r p is the constant parameter, describing the characteristic 
time of magnetic flux penetration profile; k — I/7. Now, inte- 
grating twice, the ordinary differential equation for the function 
4>(z) with the boundary conditions (45) and taking into account 
(49), we find the following explicit solution for the electromag- 
netic field distribution 



§1.5. Adiabatic approximation 

Let us now study the space and time evolution of small ther- 
mal and electromagnetic perturbations within adiabatic approx- 
imation, taking into account nonlinear flux problem. To find an 
analytical solution of equations (6) and (7) we use simple adi- 
abatic approximation, assuming that r < 1. Then eliminating 
the variable 0(x,t) by using the relationship (6) and substitut- 
ing into (7), we obtain a second-order differential equation for 
the distribution of small electromagnetic perturbation e(x, t) in 
the form 



d^e 
dz 2 



6 j P e - 

dr 



(44) 



Here, we introduced the following dimensionless variables 
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4tt j c L 2 
c 2 E c 
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1 - n 



e(z,r) = 



D = 



D 



n 2 2' 



2 2tt 
— cos —z 

T L* 



L* = 



l/ 7 



(50) 



VP' 



The obtained solution (50) describes the distribution of the 
electromagnetic field in the flux creep regime with a power-law 
current-voltage characteristics. The solution describes an blow- 
up-type instability in the superconductor sample. As easily can 
be seen that the solution remains localized within the limited 
area x < L* /2 with increasing infinitively of time. In other 
words, the growth of the solution, becomes infinite at a finite 
time t p . Typical distributions of the electric field e(z, t) deter- 
mined from analytical solution (50) is shown in Figure 1 for the 
values of parameters r p — l, n=2, j3 ~ 1 and L* ~ 0.01. 
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profile at different times for t p =1, n—2, j3 ~ 1 and L* 



0.01. 



Fig. 4c. The time evolution of the electric field profile at t p =1, 
n=2, /3 ~ 1 and L* ~ 0.01. 

Notice, that for the Kim- Anderson model the solution of the flux 
creep problem has the form 



e(z,r) = 



1 



2tt 



(51) 



Alternatively, this solution can be obtained from (50) in the lim- 
iting case, when n — > oo. 



§1.6. Nonuniform temperature profile 

It is important to notice that the profile of temperature 
T(x, t) may also significantly influence the condition of occur- 
rence of flux jumps in superconducting sample. In many cases in 
studying the dynamics of thermal instabilities of superconductors 
were usually assumed that the spatial distribution of tempera- 
ture is homogeneous in the cross section of the sample. However, 
in reality, the temperature profile of superconductor may be in- 
homogeneous along the sample as well as in its cross sectional 
plane. Such inhomogeneities can appear due to different phys- 
ical reasons. First, the vortex structure can be inhomogeneous 
due to existence of weak bonds in the sample. Second, the in- 
homogeneities may be caused by the dependence of the critical 
current density on magnetic field, the differential conductivity 
and the heat conductivity. It has been recently shown [65] that 
under the quasi-adiabatic conditions the temperature profile can 
be essentially inhomogeneous, which strongly affects the stabil- 
ity criterion of the critical state of superconductors and the first 
flux jump field. 

In this section we develop a phenomenological description of 
flux jumps in the presence of a temperature gradient and electric 
field in hard superconductors. It has been shown that thermo- 
electric effects can significantly influence the stability conditions 
in the mixed state of superconductors at very low temperatures. 
Theoretical and experimental investigations show [66-79] that in 
the mixed state during the motion of vortices inside the super- 
conductor sample may occur thermoelectric effects - temperature 
gradients. Under certain conditions the presence of a tempera- 
ture gradient in type II superconductor sample leads to critical 
state instabilities. Many investigations have been carried out 
on the thermoelectric effects in type-II superconductors, as well 
as in high-r c superconductors both experimentally and theoret- 
ically [66-79]. The discovery of high-T c superconductors stim- 
ulated intensive study of the transport phenomena caused by 
the flux-flow under an electrical current and temperature gradi- 
ent. In the presence of temperature gradient VT an additional 
thermal force appears in the sample due to the existence of a 
transport entropy of vortex lines 

F t = -sVT, (52) 
where s — s(T) is the transport entropy of a vortex line [79] 

s = s (T/T c )(l-T/T c ). (53) 
Here so is the entropy density neglecting the transport entropy of 
the flux lines. It is notice able that the existence of the transport 
entropy s is associated with the presence of low-energy electron 
states localized in the core of vortices and the value of s can be 
approximated from the local density of states [78, 79] 
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(54) 



Here e/ is the Fermi energy, fcs is the Boltzmann constant, £ is 
the coherence length. 

The critical state equation with an account the transport 
entropy has the form 



j=j(T,H,E) + SVT, 



(55) 



where S 



An additional heat released due to the transfer 



of transport entropy by a vortex lines, which is proportional to 
S. Then the equation to heat flux density q can be written as 



ST - - 
q = - K VT+—[E,B\. 



(56) 



In the quasi-stationary approximation, terms with time deriva- 
tives can be neglected in Eqs. (l)-(3). This means that the heat 
transfer from the sample surface compensates the energy dissi- 
pation arising in the viscous flow of magnetic flux in the system 
with an effective conductivity 07. In this approximation, the 
solution to equation (2) has the form 



c 



(57) 



Upon substituting the expression (55) and (57) into (6), we get 
an inhomogeneous equation for the temperature distribution 
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Here we introduced the following dimensionless variables 
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Taking into account the thermal boundary conditions (8) an ex- 
act solution to equation (58) is obtained 



Q(y) = exp 



[C 1 D v (y) + C 2 D- v - 1 (y) + eo(y)}, (59) 
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where C\ and C2 are integration constants, which are determined 
with the help of the thermal boundary conditions. In order to 
estimate the maximum temperature in the sample we present a 
solution to (58) in the form 



e(:r) = e„ 



Po- 



) 2 



near the point at which the temperature is a maximum, x = x m . 
Substituting the last relation into equation (58) we obtain an 
expression for the maximum heating due to magnetic flux jumps 



e m = 



jc+ ^^(L-Xr, 
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ckT 



(L - x n 
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(60) 



where 7 ~ 1. The calculation of the maximum heating is demon- 
strated that in the isothermal case, when the sample surface is 
cooled intensively w = woL/n 3> 1 it is easy to verify that the 
maximum heating is considerably small 6 m <§C 1 for the typical 
values of physical parameters. In the case of weak sample cooling 
id C 1 the maximum heating is about O m ~ 0, 5 -j- 2. On the 
other hand, the analysis of solution (58) show that their contri- 
bution is small in the region of x m < x < L due to the of small 
amplitude of the background electric field E = B e (x~x m )/c. On 
the other hand, near the point x = x m the temperature gradient 
dT/dx is comparatively small. Therefore, thermoelectric effects 
can significantly change the temperature profile in the critical 
state of the sample at the surface layer > x > x m only, where 
the background temperature gradients and the background elec- 
tric field are greater. Near the point of x = x m — L/2 thermo- 
electric effects are absent. The dimensionless transport entropy 
parameter p can be presented in the form 



so 
aL 



aB~ B e t K 



We B e j • (61) 

It is easily seen that p ~ I near the threshold for a flux jump, 



1/3 



when 



aBi 



I and 



B e t K 



1. The temperature rise during 



-i-Kvjc """ B e 

the flux jump depends mainly on the sweep rate of the exter- 
nal magnetic field and heat capacity and thermal conductivity 
parameters in the sample. As seen the temperature dependence 

vL 2 

p{T) is determined by the relationship t K = . At low tem- 

k(1 ) 

peratures the transport entropy increases with temperature as 

p(T) ~ S0 T~ 1/3 , (62) 

because k ~ T. This means that the influence of the thermo- 
electric effects on characteristics of the critical state stability in 
the sample can be significant under extremely low temperatures 
(T < 0, 1 K). Unfortunately, quantitative data for the transport 
entropy coefficient so are available only in limited temperature 
intervals thus the exact numerical evaluation of p proves to be 
difficult. 

Let us investigate the stability of the critical state with re- 
spect to small thermal ST and electromagnetic SE fluctuations 
with an account thermoelectric effects in the quasi-stationary 
approximation. We present solutions to Eqs. (l)-(4) in the form 



ru.t) = { ; v } <-) ( [ ) , 



E(x,t) = E(x)+ex V {f\e{ X l ) 



(63) 



where T(x) and E(x) are solutions to the unperturbed equa- 
tions obtained in the quasi-stationary approximation describing 
the background distributions of temperature and electric field 
in the sample and A is the eigenvalue of the problem to be 
determined. The instability region is determined by the con- 
dition that ReA > 0. From solution (63), one can see that 
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the characteristic time of thermal and electromagnetic pertur- 
bations tj is of the order of t K /\. Linearizing Eqs. (6)-(7) for 
small Q/T(x), e/E(x) <C 1 perturbations we obtain the following 
equations in the quasi-stationary approximation 



v\Q = k^-H + s E^ + [j + a f E] e - aEQ , 
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dx 2 
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(64) 



Eliminating the variable e between two equations in (64), we 
obtain a fourth-order differential equation for the temperature 
distribution 



we easily obtain 
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where we have used the equality 
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and boundary conditions. The right-hand side of Eq. (68) has a 
minimum at A = A c : 
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(70) 



Ar 



1 drj df 
<jfE dz dz 

ldf 



dz 



■P|i*) 



1 d 2 / 

7^ 



2 (A + 0)- 



2 f dz dz 



/0 



A 



]e = o. 



Here, we introduced the following dimensionless variables and 
parameters 
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Let us first consider the development of thermomagnetic instabil- 
ity in the adiabatic approximation, which is valid for hard super- 
conductors with t«1. In this limiting case, as seen from (63), 
the characteristic times tj of temperature and electromagnetic 
field perturbations have to satisfy the inequalities t K 3> tj 3> t m 
and At < 1 or A > 1. Using this approximation we reduce the 
last equation to a second order differential equation 
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Using the substitution 
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Ar 9 = 0. 
(65) 

(66) 
(67) 



Multiplying equation (67) by 9 and integrating the result with 
respect to y over the interval 

0<y<yi = \ f 2 (y)dy, 



equation (65) can be rewritten in the following form 
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Now we try, to obtain an integral estimation of the increment 
A, upper and lower limits of its occurrence. The behavior 
of the integrand in (70) is basically determined by the factor 
B L 

E = — — (1 — z), which is equal to zero at z = 1 (the other 

factors change more smoothly). Hence, the integrand reaches its 
maximum at z = and the upper estimate for A c is 

Xc<f\j(0),T(0),E(0)] x 1. 
It is evident that A c 1 and A c r <§C I at r <§C 1. Numerical 
evaluation gives A c « f -f- 10 2 at r = 10~ 3 . Equations (68) and 
(71) enable one to write the instability occurrence criterion in 
the form 
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Here we introduced a normalized unit vector 
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and used the inequality 
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o o 
which has been obtained with the help expansion of the function 
9(j/) in Fourier's series 



Q(y) = A m cos 



y(2m + 1)tt 



(74) 



As can bee easily seen that the shape of inequality (71) essentially 
depends on the type of coordinate dependence of the functions 
j(y) and / _1 (y). On the other hand, terms which is proportional 
y/r in (70) can be as significantly small in the limit r«l. So, 
taking into account (74) and above mentioned suggestions, the 
inequality (71) can be rewritten in the form of 
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The last inequality it may be intensify to the left-hand by using 
the relationship 
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Thus, the last stability criterion can be written in the form 
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The integral criterion for the critical state instability (76) unlike 
the analogous criterion for a homogeneous temperature profile, 
takes into account the influence of each part of the sample and 
therefore temperature gradients on the threshold for occurrence 
of flux jumps. Similar results can be found for the case of com- 
posite superconductors (t>1) when instability develops under 
a frozen-in conditions of magnetic diffusion. 

§2. Branching flux instabilities 

The dynamics of vortices in type-II superconductors exhibits a 
wide variety of instabilities of thermomagnetic origin. Nonuni- 
form magnetic flux penetration in superconductors, creating fin- 
ger and dendritic patterns, has recently attracted considerable 
interest. It is generally accepted that the nonuniform penetra- 
tion of the magnetic flux is a thermomagnetic effect due to the 
local overheating produced by the dissipative motion of vortices. 
As a consequence of the increased local temperature, the pinning 
barrier is lowered, leading to a large-scale flux invasion and to a 
final nonuniform magnetic flux distribution. Such dendritic type 
patterns driven by the flux jumping instability have been directly 
observed under a wide variety of conditions in a large number of 
superconducting samples by Duran et al. [80] and Vlasko-Vlasov 
et al. [81] in Nb films, by Johansen et al. [82], Bobyl et al. [83], 
Barkov et al. [84] in MgB2, Leiderer et al. [85] in other supercon- 
ducting materials by Bolz et al. [86, 87], Welling et al. [88] and 
Rudnev et al. [89] by means a magneto-optical imaging with a 
high spatial and temporal resolution. The existing experimental 
data [80-89] and the recently developed theoretical models [90- 
93] , suggest that the origin of these patterns is thermomagnetic 
instability of the vortex matter in the superconducting films. 

The above mentioned irregular flux avalanches and dendritic- 
like patterns has a threshold applied field and temperature, when 
the first avalanche occurs B t h, which strongly depends on both 
temperature and the sample size, as in conventional uniform flux 
jump instabilities [8-11]. The recently observed very much like 
some snow-avalanches - huge compact avalanches by Welling et 
al. [88] in Nb thin films by magneto-optical experiments also 
has a thermomagnetic origin as proposed by Aranson et al. [91]. 
The authors have determined the first flux jump field Bj as a 
function of temperature T and found that with increasing tem- 
perature more branching of the avalanches resulting in a more 
irregular flux pattern. At higher temperatures, the number of 



avalanches decreases and more flux penetration starts to dom- 
inate the behavior. However, above certain temperature these 
flux avalanches were absent. Bolz et al. [87] have studied 
the dendritic flux patterns in superconducting YBCO films by 
means a combination of magneto-optics and pulsed laser irradia- 
tion with high spatial and temporal resolution in the micrometer 
and nanosecond ranges. They found that dendrites only develop 
for certain values of the external field and temperature. Bobyl 
_et al. [83] have detected in their magneto-optical experiments 
both mesoscopic with the smallest size (50$o at B e =4 mT) and 
the macroscopic of dendritic shape with the large size (10 6 $o 
at B e =9 mT) flux jumps in MgB2 films. The results of their 
magneto-optical experiments have shown that the jump size in- 
creases with increasing of the applied field B e . Both types of 
jumps disappear above the threshold temperature (T=10 K) due 
to fast increase of the specific heat and decrease of the critical 
current. The authors believed that both types of jumps has a 
thermomagnetic origin and its existence can be explained due 
to strong demagnetization effects in thin films. Barkov et al. 
[84] believed that the dentritic structures can be observed in 
thin films with a sufficiently small thickness and in an adiabatic 
conditions, when the heat diffusivity is much smaller than the 
magnetic diffusivity. They experimentally measured a value of 
threshold field Bj in a superconducting thin MgB2 films using 
magneto-optical imaging technique which in well agreement with 
the existing theoretical data [8, 9]. It is known, that dendrite 
propagation in thin films shows velocities up to 160 km/s [86], 
i.e., these velocities are much higher than the speed of sound. 
Wertheimer and Gilchrist discovered a well defined pattern of 
flux dendrites with a propagation velocity v in the interval be- 
tween 5 m/s and 100 m/s [94]. The dendrites velocity depends 
on the disks thickness, for smaller d a higher velocity was found. 
Bolz et al. [86] have found that for thin superconducting films 
the flux jump instability can give rise to dendritic magnetic flux 
avalanches propagating with tip velocities as high as 50 km/s. 
The dendritic type of flux penetration into type-II superconduc- 
tor slab was proposed by Aronson and his co-workers [91] on 
the basis of numerical simulations of coupled equations for the 
magnetic induction and temperature in the limit of weak Joule 
heating in the sample. A formation of such dendritic flux front 
patterns can be observed under condition that the magnetic flux 
diffusion is much faster than the heat diffusion. They proposed 
that each vortex microavalanche results in a partial flux penetra- 
tion process which is accompanied by local Joule dissipation. It 
has been pointed also that the dendritic microavalanches may be 
also initiated by macroscopic surface defects, which can trigger 
a global flux jump instabilities in high-T c superconductors. 

§2.1. Flux flow 

Let us consider the problem of occurrence branching instabilities 
in the flux flow regime. For this purpose we formulate a ba- 
sic equation governing the dynamics of magnetic field induction 
and temperature perturbation in a semi- infinitive superconduct- 
ing sample x > 0. The distribution of the magnetic flux density 
B and transport current density j inside a superconductor is 
given by a solution of the equation 



dB 
dx 



VoJc 



(77) 



The electric field E(x, t) is generated inside the sample according 
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(78) 



to Faraday's law 

dE _ dB 
dx dt 

In the flux flow regime the electric field E(x, t) induced by the 
moving vortices is related with the local current density j(x,t) 
by the nonlinear Ohm's law 



E = p(j-j c ). 



(79) 



Let us suppose that in the flux flow regime the differential re- 
sistivity is approximately constant and independent on magnetic 
field, i.e. p = p/=const [41]. The spatial and temporal evolution 
of thermal T(x, t) and magnetic field B(x, t) perturbations are 
described by the thermal diffusion equation coupled to Maxwell's 
equations 

de rf 2 e +2 ^ + a 
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dz 2 



dz 



(80) 
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where we have introduced the dimensionless variables and pa- 
rameters 
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We present the small thermal and electromagnetic perturbations 
in the form 



Q(z, t) = 6 exp [jt/to + ik£ + iq(] , 
e(z, t) = e exp [it /t + ik£ + iqQ , 



(81) 



where 7 represents the eigenvalue of the problem to be deter- 
mined and having wave numbers k,q = 2nd/L in the £ = x/L 
and C = y/L directions, respectively. Solving the system equa- 
tions (80) by using the solution (81) we obtain the following 
equation for the eigenvalues 

7 2 + [(l + r)(fc 2 + q 2 )-l] 1 +[r(k 2 + q 2 )-l](fc 2 + q 2 )-2(k + q) 2 = 

where r is the ratio between the characteristic time of magnetic 
flux diffusion and the characteristic time of heat flux diffusion 
[4]. The instability criterion for flux jumping is determined by a 
positive values of the growth rate Re 7>0. Let us first consider a 
simplest case, when r=0. Then, from the last relation we obtain 



7 2 + [{k 2 + q 2 ) - 1] 7 - (k 2 + q 2 ) - 2(k + qf 



0. 



(82) 



First, we notice that for perturbation uniform in x-direction k~0 
the growth rate is positive Re7 > for wave number q < 1. In 
such case the small perturbations grow with the maximal pos- 
sible rate Re7 = 1. For the case q > 1, the growth rate 7 
is negative, consequently, the small perturbations always decay. 
Similarly, for the perturbations uniform in y-direction q=0, the 
growth rate is positive Re7 only, if wave number is satisfied k<l. 
As the wave number q approaches infinity q — > 00, the growth 
rate approaches 7=1 for k — and small perturbations grow 
with time. It can be shown that, for the critical values of the 
wave number k = k c , the growth rate is zero Re7 = 0. This 
instability occurs at k < k c> and its criterion can be written as 

k < k - — 



The growth rate 7 dependence on the wave number q for differ- 
ent values of k are demonstrated in Figs. (5a-5d) at r = 0.1. For 
high enough values of wave number k the system is stable. As 
the wave number k decreases, the growth rate 7 increases. The 
branching instability will gradually appears for relatively small 
values of the wave number k=0.5. 




2 3 4 

Wavenumber Q 

Curve 1 k=0.1 

Curve 2 k=0.3 

-■ Curve 3 k=0.5 



Fig. 5a. The dependence of growth rate of the wave number 
k=0.5, 0.3 and 0.1 for r = 0.1. 



Thus, on the basis of a linear analysis of a set of differ- 
ential equations describing small perturbations of temperature 
and electromagnetic field we found that under some conditions 
a branching instability may occur in the sample. 
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Fig.5b-5d. The dependence of the growth rate on the wave 
number for k=0.5, 0.3 and 0.1 for r = 0.1. 



§2.2. Flux creep 

Let us now study qualitatively, the problem of occurrence 
of the branching instabilities within the framework of flux creep 
problem with a nonlinear voltage-current characteristics. We as- 
sume that an applied field parallel to the surface of the sample. 
As we have showed above, for the flux creep problem the differ- 
ential conductivity a is determined by the following expression 



dE 



3c 

nE b 



(83) 



Taking into account this flux creep relation, we write the system 
of differential equations, governing the small perturbations of the 
temperature and electromagnetic field in the form 
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(85) 



Next, we shall present a solution of these differential equations 
for Q(x,t) and electromagnetic field e(x,t) perturbations in the 



following form 

6(3!, t) = To(x) + (T c - To)9 exp [jt/to + iqz] . 

e(x,t) = E b (x) + _Ei,eexp [jt/to + iqz] . 



(86) 



where To (a:) and Eb(x) are solutions to the unperturbed equa- 
tions obtained in the quasi-stationary approximation describing 
the background distributions of temperature and electric field in 
the sample. From solutions (86), one can see that the charac- 
teristic time of thermal and electromagnetic perturbations t is 
of the order of to/7- Here, we have introduced the following 
dimensionless parameters and variables 



7T t 
H' 



to = ^— , z = j, 1= — -, q = 

As we mentioned above, the background temperature Tq(x) is 
practically uniform over the cross-section of the sample and un- 
der this approximation we ignore its coordinate dependence. It 
turns out that these simplifications have no qualitative influence 
on the results but make it possible to perform analytical calcu- 
lations completely. 

Substituting the last expression (86) into the system equa- 
tions (84), (85) one can get the following linearized system equa- 
tions for O and e 

rq 2 Q + 1 Q+ -@-2 (-'] e = 0, 

W (87) 

q 2 e + 7 [e - 



n@] = 0, 

Solving the last system equations we obtain the following disper- 
sion relation for the eigenvalue problem 
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The instability of the flux front is defined by the positive value of 
the rate increase Re 7>0. It can be seen that there is a critical 
cutoff wave number, 

4c = ^. (89) 

below which the system is always unstable at n=l. This insta- 
bility appears first at q=0. In this case the small perturbations 
grow with the maximal possible rate 7 = 1. The growth rate 
dependencies on the wave number for different values of r are 
demonstrated in Figs. 6a-6c at n = 1. For high enough values 
of t the system is stable. As the r decreases, the growth rate 
7 increases. The branching instability will gradually appears for 
relatively small values of r=0.5. 

According to inequality (89), the branching instability occurs 
at the threshold electric field E = E c , which can be written as 



E r = 



7T k(Tc — T ) 



nj c L 2 



Taking into account that, the penetration depth L, the threshold 
field can be written at B = B t h as 



Bth = 



k(T c — T )j c 
nE b 



The threshold field for branching instability, as can be seen from 
the last expression is highly sensitive to the critical current den- 
sity and the shape of the background electric field E b , generated 
by the varying of magnetic field. The threshold field Bj decreases 
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monotonously with increasing the background electric field Eh. 
If we assume that the background electric field E\, generated by 
a varying of magnetic field as Eb — B e , for the considered simple 
geometry, then we can easily obtain the expression for the sweep 
rate dependence of the threshold field of branching instability. 



Let us assume that the thermal diffusion is slower than the 
magnetic diffusion t < 1. In this limiting case, the instability 
criterion is determined as q = q c = 1 for n=l, so the threshold 
field can be presented as 



B t , 



Mo 



-{T c -To). 



Wavenumber q 



This is a well known adiabatic criteria [2], which assumes that 
the heat transport from the sample surface to the environment 
can be neglected. 
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Fig.6a-6c. The dependence of the growth rate on the wave 
number for r = 0.5, 0.05, 0.08. 

Fingering-like patterns were found theoretically recently by 
Rakhmanov et al. [90] on the linear analysis of a set of thermal 
and Maxwell equations for small temperature and electric field 
perturbations within the framework of adiabatic approximation. 
They calculated the stability criterion and estimated the build- 
up time and finger width. The fingering instability occurs at 
larger electric fields E > E c , only. The authors believed that 
such fingering instability may develop to dendritic flux penetra- 
tion into the sample. By solving the Maxwell and the thermal 
diffusion equations, it was shown that for small fields there are 
no solutions for perturbations growing in time, implying a sta- 
ble situation. As the field increases the distribution can become 
unstable, with a fastest growing perturbation having a non-zero 
wave vector along the film edge. This means an instability will 
develop in the form of narrow fingers perpendicular to the edge - 
a scenario closely resembling the observed dendritic flux behav- 
ior. 

§3. Dynamically driven flux instabilities 

Recently, a many small in magnitude - mesoscopic in nature 
flux jumps has been observed in magnetization measurements of 
conventional and high-T c superconductors at very low temper- 
atures. It is interesting that the appearance of such type flux 
jumps can not be described in the framework of the traditional 
model of magnetothermal instabilities [1-9]. The main features 
of observed flux avalanches are the following: the first flux jump 
occurs at a values of the applied field much higher than full pen- 
etration field B p ; the observed avalanches of almost the same 
amplitude look quite regular, they are essentially magnetic his- 
tory sensitive; the avalanches exist only above threshold field of 
several tesla in magnitude as opposed to the low-field limit; its 
onset is independent on the sweep rate of the external magnetic 
field, a distinct contrast to flux jumps originating in a thermal 
instability in type-II superconductors; the avalanche size distri- 
bution is strongly peaked with a characteristic size, clearly dis- 
tinct from the broad power-laws in the self-organized criticality 
concept; the jumps occur at significantly different fields on the 
increasing and decreasing branches of the hysteresis loop. Thus, 
it can be concluded that the onset of macroscopic flux jumps 



16 



in these experiments may be driven by a dynamic instability, 
which may be explained by means of a concept of self-organized 
criticality. 

For example, Seidler et al. [95] argued for a dynamical origin 
of small vortex avalanches observed in their magnetization ex- 
periments and proposed an explanation based on self-organized 
criticality (SOC) [96, 97]. A new type of macroscopic flux jump 
with a narrow size distribution was observed by Zieve et al. [98] 
in untwined single crystals at very low temperatures (well be- 
low T ~1K), one not triggered by thermal instabilities. Inter- 
estingly, the avalanches exist only above threshold field and its 
onset is independent on the sweep rate. These avalanches occur 
at significantly different fields on the increasing and decreasing 
branches of the hysteresis loop. The authors concluded that the 
flux jumps in their experiments may be driven by a dynamic 
instability, which may be explained by means of a concept of 
self-organized criticality. Based on MD simulations Olson et al. 
[99] have showed that the dynamical instabilities are triggered 
when the external magnetic field is increased slightly, and are 
thus driven by a flux gradient rather than by thermal effects. 
The existence of SOC type of behavior flux avalanches recently 
has been experimentally observed and numerically adopted by 
many researchers [100-109]. Altshuler et al. [109] recently have 
analyzed in detail an experimental results on the flux avalanches 
which may be qualitatively understood by the concept of a self- 
organized criticality. 

Huge "catastrophic", flux-jump-like avalanches associated 
with sudden movement of many vortices were observed by Be- 
hina et al. [110] in a superconducting Nb thick film during a slow 
sweep of external magnetic field (B e ~1.1 Oe/s) by means a Hall 
probes. The obtained in a series of measurements avalanche size 
statistics indicated that the size distribution of these avalanches 
presents a power-law behavior only a limited range - in the small 
event region. The authors believed that in contrast, at low tem- 
peratures and fields, huge avalanches which may be related to 
thermal instability of the Bean state dominate the dynamic re- 
sponse of the sample. 

Irregular and non-periodic flux jumps were observed by Ma- 
jumdar et al. [Ill] in the low field magnetization measurements 
at different field sweep rates in the heavy fermion superconduc- 
tor samples. The large in magnitude flux jumps were observed 
at the highest sweep rates. The magnitude of these jumps de- 
creases slowly with decreasing of temperature below 1.7 K and 
the jumps are completely absent at 1.5 K. They argued that the 
observed flux jumps are due to local flux entry through the sur- 
face or geometrical barriers. The asymmetry of the flux jumps 
with respect to increasing and decreasing part of the magnetiza- 
tion loop clearly indicate the existence of barriers at surface of 
the sample. 

Flux jumps that differ qualitatively from well-known magne- 
tothermal instabilities have been observed by Milncr [112] and 
Gerber and Milner [113] in thin single crystals of high-T c super- 
conductors in the fields up to 17 T at the temperatures down to 
0.3 K, using magnetometric, calorimetric and induction pick-up 
techniques. The flux jumps start above full penetration field B p , 
demonstrating clear periodicity in the slowly changing magnetic 
field. The magnetization hysteresis loop is non-homogeneous, 
consisting of high steps, which reflect the global redistribution 
of the magnetic flux in the sample. It was found that the fre- 
quency of flux jumps is strongly temperature dependent. Milner 



proposed a few variants to explain about the origin of the low- 
temperature/high-field magnetization jumps observed in their 
experiments. Using a Hall probe sensor Terentiev et al. [114] 
have observed a quasi-periodic flux instabilities below a certain 
temperature (T ~3 K) in a superconducting Nb thin films by 
a square lattice of Ni dots. The magnetization measurements 
indicated that the Ni dot lattice exerts a crucial influence on 
the appearance and nature of the instabilities. On the analysis 
of magnetization curves the authors have found that the quasi- 
periodic instabilities as unexpected low-temperature matching 
anomalies, most probably initiated near the film edge, where 
flux density is much lower than in the film interior. 

§4. Flux avalanches 

Recently, Chabanenko et all. [26] have reported an interest- 
ing phenomenon in their experiments - convergent oscillations 
of the magnetic flux arising from flux jump avalanches. The au- 
thors argued that the observed oscillations due to flux avalanches 
can be interpreted as a result of the existence of a definite value 
of the effective vortex mass [27]. Thus, it is necessary to take 
into account collective modes, i.e., the inertial properties of the 
vortices in studying the dynamics of the flux avalanches. Prior 
to the jump, the mixed state of superconductors is characterized 
by nonuniformly distributed magnetic induction localized near 
the surface. As a result of the avalanche, the flux rushes from 
either sides of the sample towards the center [26]. Two fronts 
of the penetrating flux collide in the center of the sample and, 
owing to the existing vortex mass, give rise to the local surplus 
density of the magnetic flux that exceeds the value of the exter- 
nal magnetic field. The repulsion force in the vortex structure at 
the center of the sample that have resulted from its compaction, 
initiates the wave of the vortex density of the inverse direction 
of propagation. Upon reaching the surface, this wave is reflected 
from it. This results in the oscillations in the vortex system [26]. 
The limitation of the number of oscillations observed is caused 
by the existence of damping. One succeeds in observing the os- 
cillation of the vortex density only owing to a strong compression 
of the vortex structure as a result of the giant avalanche-flux [17, 
24-26]. 

In this section, we study the dynamics of the magnetic flux 
avalanches, which take account inertial properties of the vortex 
matter. 

In the flux flow regime the electric field E(r, t) induced by 
the moving vortices is described 

E = vB. (90) 

To obtain quantitative estimates, we use a classical equation of 
motion of a vortex, which it can derived by integrating over the 
microscopic degrees of freedom, leaving only macroscopic forces 
[32]. Thus, an equation of the vortex motion under the action of 
the Lorentz, pinning, and viscosity forces can be presented as 

m^+vV + F L +F p = 0. (91) 
at 

Here m is the vortex mass per unit length, Fl = -j$o is the 

c 

Lorentz force, F„ = -j c $ , V is the flux flow viscosity coefficient 

c 

[32]. For simplicity we have neglected the Magnus force, assum- 
ing that it is much smaller then the viscous force (for example, 
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for Nb see, [26]). In the absence of external currents and fields, 
the Lorentz force results from currents associated with vortices 
trapped in the sample. 

In combining the relation (90) with Maxwell's equation (2), 
we obtain a nonlinear diffusion equation for the magnetic flux 
induction B(r, t) in the following form 



dB _ 
dt 

dV 

m— — + r/V — 



V(uB) 
1 



$o(j ~je), 



(92) 



(93) 



dt c 

The temperature distribution in superconductor is governed by 
the heat conduction diffusion equation 

v(T)^ = V[ K (T)VT\+jE, (94) 

Let us present a solution of equations (90-94) in the form 
T = T + e(x,t), 

B = B e + b(x,t), (95) 
V = V +v(x,t) 

where To(x), B e (x) and Vo(x) are solutions to the unperturbed 
equations, which can be obtained within a quasi-stationary ap- 
proximation. Substituting the above solution (95) into equations 
(92-94) we obtain the following system of differential equations 
[26] 



, = 2v 
dt 

dv 
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dt \dx J \dx 
where we have introduced the dimensionless variables 
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Assuming that the small thermal and magnetic perturbations has 
form Q(x,t),b(x,t),v(x,t) ~ exp[7t], where 7 is the eigenvalue 
of the problem to be determined, we obtained from the system 
Eqs. (96)- (98) the following dispersion relations to determine the 
eigenvalue problem 



(7 + /3) 



d 2 b 
dx 2 



[(7 + 0)n - 2/3] 



db 
dx 



[(M + l)7 2 + 



+ [(p - 1)0 -j»-1] 7 -0*- 1)06 = (99) 

An analysis of the dispersion relation shows that, the grows rate 
is positive Re 7>0, if fi > fj. c = 2 and any small perturbations 
will grow with time. For the case when y, < fi c , the growth rate 
is negative and the small perturbations will decay. At the critical 
value of n = fi c , the growth rate is zero 7=0. For the specific 



case, where y, = 1 the growth rate is determined by a stability 
parameter p. Thus, the stability criterion can be written as 



P>1. 



(100) 



For the case, where thermal effects is negligible (j3 = 1) we may 
obtain the following dispersion relation 



d?b 

dx 2 



+ (7- !)(/* + 1)& = 0. (101) 

dx 



Seeking for b ~ exp(ifca;) in dispersion relation, the growth rate 
7 dependence can be obtained as a functions of wave number k. 
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Fig. 7a. The dependence of the growth rate on the wave number 
for n = 0.1,0.5,0.8. 

The stability of the system depends on the growth rate, 7, given 
in (101). We analyze the growth rate of small perturbations as 
a function of wave number k. When k < k c — n the growth 
rate is positive and any small perturbations will grow with time. 
For wave number k > k c , the growth rate 7 is negative. Conse- 
quently, the small perturbations always decay. It can be shown 
that, for wave number k = k c the growth rate is zero 7 = 0. As 
the wave number approaches zero k — > or infinity k — ► 00 
the growth rate approaches 7 = 1 and small perturbations grow 
with time. As the wave number approaches unity k = 1 the 
growth rate is determined by the value of ji 



7 ; 



2fi 



For fi = the growth rate is zero 7 = 0. For fj, = 1 the growth 
rate is unity 7 = 1. Since the growth rate is zero at the critical 
wave number and approaches to unity in the limit of zero wave 
number, there must exist a wave number in between that max- 
imizes the growth rate. Figs. (7a-7d) show the growing rate, 
7, as a function of the wave number k, for different values of 
the parameter /i. As the value of fi increases, the corresponding 
growth rate increases. 
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Fig.7b-7d. The dependence of the growth rate on the wave 
number for fi — 0.1, 0.5, 0.8. 



Thus, we show that at under some conditions flux jump 
avalanches may occur in superconductor sample, which take into 
account inertial properties of the vortices. 



§5. Second magnetization peak 

The second-peak effect in the magnetization curve is one of 
the most remarkable features in the vortex state of both low 
and high-r c superconductors. The second magnetization peak 
(or "fishtail") as a function of applied field was observed in a 
low-temperature superconductor materials [115-118] and high- 
temperature superconductors [119-130]. The second peak is ex- 
pected to give an important clue in the understanding of the 
complex vortex-matter phase diagram in low and high-T c super- 
conductors. In addition, different interpretations for the second 
peak have been discussed in literature, however, the mechanism 
of the second peak remains unclear. These range from the en- 
hancement in pinning due to matching effects in oxygen defi- 
cient structures [119], the collective creep phenomenon [120, 130], 
matching effect [131, 132], the dimensional crossover [133], the 
surface barrier effect [134], the thermomagnetic instability [135- 
139], elastic to plastic creep [122-124], crossover in the pinning 
regimes, e.g. from single-vortex pinning to a pinning of vortex 
bundles [121]. Many theoretical and experimental publications 
dealing with the second magnetization peak effect have been re- 
cently published, the interpretation of this effect is still rather 
controversial. No single mechanism exist until day, which may 
explain the nature of these effects that were observed in a large 
amount of superconducting samples. 



§5.1. Thermomagnetic instability 

For low-T c superconductors (Nb), the second magnetization 
peak occurs at lower temperature region during the magnetic 
flux penetration, but at higher values magnetic fields close to the 
upper critical field H C 2- Resent experimental investigations indi- 
cated that second magnetization peak is probably originated by 
a thermomagnetic instabilities that occur in the low-temperature 
region, far below H C 2 [135-139]. In particular, such flux-jump in- 
stabilities in Nb thin films have been studied by Esquinazi et al. 
[136] on the basis of isothermal and global magnetization mea- 
surements as a function of applied field with a SQUID and /i-Hall 
sensors, respectively. Large jumps as well as pronounced irreg- 
ularities of large scale gradients in the vortex density has been 
observed at low temperature region (T ~ 3 K) and at low ap- 
plied fields (H ~ 100 Oe) in the magnetization M(H) curve. The 
existence of the second magnetization peak in a thick Nb film, 
which was clearly produced by magneto-thermal instabilities re- 
cently have been reported by Kopelevich et al. [137, 139] using a 
SQUID and Hall array measurements. The magnetization peak 
appears at fields higher than the first maximum observed usu- 
ally in a zero-field-cooled state, but much lower than the upper 
critical magnetic field Hc2. The flux jumps were clearly ob- 
served up to the second magnetization peak at B sp « 1300 Oe. 
The first flux jump field Bfj was discussed and numerically es- 
timated in this work. This estimation shows that the first peak 
Hf p ~ 100 Oe observed in the magnetization curves is of the or- 
der of the first flux jump field Bfj. Additionally, the temperature 
and sample size dependencies on the magnetization jumps were 
investigated. It has been shown that the maximum dissipated 
magnetic energy during the flux jump process - the maximum 
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temperature at the first flux jump may increase substantially at 
B > Bfj due to the decrease of the specific heat and the in- 
crease in the critical current density in the sample. Thus, under 
certain experimental conditions, which depend on the thermal 
and magnetic diffusivities and the thermal coupling of the sam- 
ple with its environment, this temperature change can produce 
flux jumps instabilities. It was also shown [136] that the field 
of the second peak increases with the thickness and/or width of 
the film as expected for measurements in transverse geometry. 
The jumps become smaller in magnitude or even vanish with a 
reduction of the film thickness. The authors have pointed out 
that careful characterization of the second peak in Nb films is 
important and may contribute to understand the anomalous be- 
havior of the magnetization curve M(H) measured in high-T c 
as well as conventional superconductors which in some cases is 
interpreted in terms of a vortex lattice phase transition [136]. 
Similar magnetization measurements on the second magnetiza- 
tion peak were performed recently Stamopoulos et al. [140] at 
different magnetic field sweep rates in superconductor thin films. 
Thermomagnetic flux jump instability was observed at temper- 
ature T—7.2 K. The magnitude of the first flux jump field Bfj 
coincides with the first peak B p j at a characteristic temperature 
T = 7.2 K. The first flux jump field remains constant Bfj Oe at 
low temperature region T < 6.4 K. The dependence of the first 
flux jump field on sweep rate of the magnetic field has been stud- 
ied, also. Surprisingly that the observed first flux jump field is 
not inversely proportional to the sweep rate of the applied field, 
which in contrast to theoretical concepts for the flux jumping 
and other experimental dates [8, 15, 136, 17]. 

Let us qualitatively estimate the temperature dependence of 
the first flux jump field _B/j(T), which occurs at low tempera- 
ture part of the magnetization curve. According to conclusion of 
Ref. [136] the first magnetization peak Bf p (T), which is nearly 
independent of the sweep rate is of the order of the field at which 
the first jump Bfj(T) occurs. It is reasonable then that the tem- 
perature dependence of the first flux jump field Bfj(T) can be 
fit according to the adiabatic critical state theory 

In order to compare the last formula with experimentally de- 
termined first flux jump field, one must know the temperature 
dependencies of the critical current density and the specific heat 
of the sample. As has been mentioned in literature [136] a quan- 
titative estimation of ic(T) from the available experimental and 
theoretical models is difficult because of the uncertainty in the 
values of the critical current j c at a given field and temperature. 
To determine j c {T), different approaches have been taken. Em- 
pirically, the critical current density J C (T) can be presented in 
the form 

j c = jo[i-n m . (102) 

where 1 < m < 2. The different exponents n=l and 2 refer to 
the most common cases discussed in the literature, where the 
critical current exhibits a linear and a quadratic dependence on 
t = T/T c . 

There is experimental evidence [41, 33], which indicates that 
the temperature dependence of the critical current density ap- 
proximately linear at low temperatures (n=m=l). On the other 



hand, at temperatures well below T c the magnitude of the spe- 
cific heat decreases rapidly with decreasing temperature, so we 
can estimate the specific heat as v ~ (T /T C ) A . Using the above 
values of temperature dependencies of the specific heat and the 
critical current density we obtain the following stability criterion 
[1, 2] 

Bfj = T^\/ n 3 T 3 (T c — T). 

Combining all these equations, one obtains a theoretical Bfj(T) 
curve as shown in Fig. 3. Thus, it can be concluded that the 
temperature dependence of the first penetration field B/ P (T) is 
consistent with the temperature dependence of the first flux jump 
field B/p(T), as demonstrated in figure 8. 
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Fig.8. The temperature dependence of the first flux jump field 

The dome-like profile of the temperature dependence of the first 
flux jump field at low fields and temperatures have been observed 
in many experiments [136, 143]. Thus, our analysis shows that 
the temperature dependence of the first flux jump field Bfj(T) 
is consistent with the model which ascribes the first peak (on- 
set field) of the second magnetization peak. This scenario is 
consistent with the observation of a crossover in the field pro- 
files across the sample, from profiles characteristic of geometrical 
barrier with weak pinning at fields below the Bf p , to Bean-like 
profiles at fields above the first peak. The results of experi- 
ments, obtained by Esquinazi [136] have shown that near the 
first flux jump field the vortices penetrate the sample, forming 
a droplet with a dome shape far away from the Nb film edges, 
i.e., near the center or at the center of the sample, depending 
on the temperature. A further increase of the external field does 
not change homogeneously the flux profile inside the sample, but 
part of the vortices remain pinned. The domelike flux profiles at 
low enough fields have been also observed in high-temperature 
platelet-shaped samples (see, for example [143] due to the ex- 
istence of a geometrical or edge barrier given by the shape of 
the sample. It has been claimed [136, 138] that the observed 
domelike shape of the field profile cannot be due to edge barrier. 
The observed strong temperature dependence of the flux profile 
at the first flux jump field Bfj(T), speaks for the influence of a 
thermomagnetic instability [136]. After the first flux jump, the 
increase of field moves vortices from the edges and a Bean-like 
profile starts to develop. The higher the applied field, the lower 
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the critical current density and therefore the larger is the pene- 
tration of the Bean-like profile into the sample up to the second 
magnetization peak field B sp . 

The results of [136] indicated that the temperature depen- 
dence of the second peak field B sp (T) has shown a similar tem- 
perature dependence as the upper critical field B C 2(T). So, at 
high enough fields and temperatures the temperature dependence 
of second magnetization peak field can be well fitted by 

B sp (T) = B (1 - t 2 ) . 



where Bo is a constant parameter. The temperature dependence 
of the field B sp (T) is shown in Fig. 9. 
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Fig. 9. The temperature dependence of the second 
magnetization peak. 



As can be seen it reflects, roughly the temperature dependence 
of the upper critical field. At higher temperatures and fields 
(near the upper critical field), where the second magnetization 
peak appears the temperature dependence of the critical current 
density can be presented as 



Jc(T)=j c (0) {l-t 2 ) 2 



commonly accepted in literature [33] (n=m=2). Assuming that 
the thermal conductivity is a linear function of temperature, we 
can easily obtain an expression for the temperature dependence 
of the field Bj (T) 



Bj (T) * Bio y*CF 



Fig. 10. The temperature dependence of the threshold field 
B 3 (T). 

It should be noted that the stability criterion (43) strongly de- 
pends on the sample geometry. For the case of infinite slab geom- 
etry or long cylindrical sample with an external magnetic field 
applied parallel to its surface or axis the role of dimension is 
played by the sample diameter d — 21. Then the critical dimen- 
sion of the sample is determined from the relation 

47T 

Bj = — j c d c , 
c 

Then for the sample with a dimension d smaller than the criti- 
cal one d < d c , no flux jumps occur at any temperature for any 
external magnetic field. Thus by reducing the sample dimension 
it is possible to avoid flux jumps. In the case of disc geometry 
the stability criterion can be modified as 



13 = 



irjcrd 
v(T c -T o y 



(103) 



where r is the radius and d is the diameter of the disc. In this 
case the experimental value of stability parameter for Nb disc 
samples approximately is equal to /3=1.2 for the typical values of 
parameters j c , v and r, d at helium temperature. For thin films 
in an external magnetic field perpendicular to their surface, the 
role of the critical dimension from the point of view of occurrence 
of flux jumps is played by their thickness. The result obtained for 
bulk sample may be modified for thin samples in a perpendicular 
field multiplying by a numerical factor. The dependence of the 
maximum in the width of the magnetization hysteresis loop and 
transition field Bj - on the sample geometry has been reported 
for Nb films and single crystals [137]. It was found that the 
second magnetization peak vanishes in both superconductors for 
small enough dimensions of the sample. According to flux jump 
theory the sample dimension is 



d = 



-B , 



where Bio is the free parameter. 



47T/ C 

This model predicts the vanishing of the second magnetization 
peak with decreasing the sample size. The estimation of the 
critical thickness of the sample from the measured values of pa- 
rameters and Bj show that the necessary requirement for the 
avoiding flux jumping the sample thickness has to be lower than 
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its critical value, of the order of 100 \im [136]. In agreement 
with this model, the second peak vanishes in thin film super- 
conductors when the lateral sample size becomes d less than 
d < dc ~ WOfim. For the transverse geometry of a rectangular 
thin film the effective size d of the sample above which flux jump 
occurs is [137] 



where d « l,w; I is the length, w is the width of the film. For 
the numerical magnitudes of a typical parameters, the critical 
seise d of the sample for the transpierce geometry should be less 
< 500/im [137]. Hence, the authors argued that the existence of 
the second magnetization peak in thick Nb film can be should be 
clearly produced by thermomagnetic flux jump instabilities [8, 
9]. The authors [136] believed a such sample size dependence of 
the threshold field as well as the observed dependence of the am- 
plitude and number of flux jumps on the field sweeping rate B e 
provide the experimental evidence for the thermomagnetic origin 
of the second magnetization peak in Nb samples. Let us estimate 
the temperature dependence of the effective sample size d c for a 
typical experimental values of parameters for Nb. In Fig. 11 we 
show the temperature dependence of the sample size d. Roughly, 
it increases linearly with temperature. This result was obtained 
from the relation for the stability criterion (43) and it can be 
considered as empirical that we assumed a linear dependence of 
the thermal conductivity k = kq(T/T c ), where kq is the con- 
stant parameter, as suggested by low-temperature data for Nb 
and linear temperature dependence of the critical current density 
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Fig. 11. The temperature dependence of the effective sample 

size. 

Nabialek et al. [17] have shown that due to the existence of 
a critical dimension the flux jump instability in high-T c super- 
conductors was observed only in relatively large single crystals. 
No flux jumps has been observed in ceramic materials, because 
the critical dimension in these materials is limited by the grain 
size, which is typically very small. Thus, the nature of origin and 
further development of the flux jump instabilities may be differ- 
ent in thin films and bulk samples [93]. In bulk samples flux 



jump can be triggered when the gradient of magnetic flux profile 
inside the sample exceeds critical value Bj. In each jump many 
vortices are involved in the avalanche process, which normally 
expands to a large part of the sample volume. On the other 
hand, in bulk samples, there are many large grain boundaries, 
which act as strong pinning centers. The gradient of the flux pro- 
file near these boundaries can be broken at a certain limit. Once 
a blast occur, the thermal energy induced by drastic flux mo- 
tion cannot easily diffuse out and be carried by the environment. 
Therefore this self-heating will lead to an increase of tempera- 
ture in the region in which the flux instability occurs, leading to a 
large jump on the magnetization. One can understand from this 
picture that the number of flux jumps cannot be large in bulk 
samples. Such instabilities in the bulk superconductor samples 
can be qualitatively understood based on the adiabatic theory 
[1]. In thin films, the situation is completely different and the 
observed small and irregular avalanches cannot be explained by 
the adiabatic theory. There are generally more defects or pinning 
centers in thin films may be observed a small avalanches. Due 
to the high density of small defects formed during the prepara- 
tion process there are many places for the avalanche to occur in 
thin films. Thermal diffusion is much easier in thin films due to 
their very small thickness and large surface area exposed to the 
environment. Therefore in thin films each avalanche is small in 
magnitude but the number of avalanches can be huge. This pic- 
ture may give an explanation to many small vortex avalanches 
observed in thin films. Since maintaining good thermal contact 
throughout a sample is more difficult for bulk samples than for 
thin films, thermally driven avalanches should occur more easily 
as sample thickness increases [93]. 

We have studied the second magnetization peak in the frame- 
work of a thermomagnetic model and demonstrated by a qualita- 
tive agreement of the theory with experimental results [135-139]. 
In other words, our theoretical analysis on flux jump instabil- 
ity qualitatively reproduces experimentally observed in [135-139] 
second peak features: the sample size dependence of the thresh- 
old field, the temperature and sweep rate dependencies of the 
first flux jump field and etc,. 

Magnetic flux jump instabilities have also been observed in 
high temperature superconducting materials [2, 18, 62, 19, 62, 
17]. In the following we shall discuss some recent existing theo- 
retical models, which may give a qualitative description on the 
nature of the second magnetization peak and its onset for su- 
perconducting materials. Both the onset field B on and the sec- 
ond peak field B 3p exhibit strong temperature dependence up 
to the close vicinity of critical temperature T c . To understand 
the origin and nature of the onset of second peak it is very use- 
ful to study the temperature dependence of the field B on . Let 
us study the temperature dependence of the onset field within 
the framework above mentioned theoretical models. In the last 
decade, many different by each other hypotheses have been pro- 
posed to explain the nature of the second peak effect. For ex- 
ample, Krusin-Elbaum et al. [121] ascribe the occurrence of the 
peak effect to the crossover between two collective vortex creep 
regimes; from single to collective vortex creep regime. The in- 
dependence of the magnetic moment from the applied field for 
low temperature isothermal loops is interpreted by them as the 
single vortex creep regime [121]. The authors believed that the 
peak effect is not likely to be observed on very short time scales 
due to the fact that the relaxation is much faster in the single 
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vortex creep regime than in the collective vortex creep regime. 
Yeshurun et al. [130] also assigned dynamical characteristics to 
the appearance of the peak effect. In other studies, the influ- 
ence of oxygen content on the peak behavior is discussed [142]. 
Zhukov et al. [142] found that the peak effect was present at 
high temperatures for all oxygen contents, however the peak's 
height takes the smallest value for the state closest to stoichio- 
metric status, which means that oxygen deficiency is necessary 
for achieving greater pinning and therefore higher critical current 
density. In some cases a correlation between the peak effect and 
the oxygen concentration is found to be the source of peak effect. 
Daumeling et al. [119] argued that a field-induced pinning en- 
hancement in oxygen-deficient regions may be the source of the 
second magnetization peak in superconductor sample. The au- 
thors proposed that oxygen-deficient areas had become normal as 
the field was increased, as a result of their lower T and H. These 
areas thus become new pinning centers. Therefore any change in 
oxygen content or distribution should lead to suppression or ap- 
pearance of the peak effect. In a quite different model, the peak 
effect is considered as a result of the matching effect [132] , where 
matching of the vortex lattice with the defect structure leads 
to pinning enhancement. However, the temperature-dependent 
nature of the peak strongly rules out the matching effect as the 
origin of the peak effect. 



§5.2. Collective creep 

One of simple situation to introduce the concept of weak 
collective pinning is an isolated vortex in an isotropic supercon- 
ductor subject to pinning by randomly distributed point defects 
[141]. We shall study the problem within a continuum elastic 
description of the flux-line lattice, with a potential energy given 
by the combination of the elastic energy, the pinning energy and 
the action of the Lorentz force 



2 r U 

C66« L + d — 



1/2 



(104) 



where u and L are the transverse and longitudinal sizes of vortex 
distortion respectively, 7 = f p ni(, is the disorder parameter; f p is 
the elementary pinning force, m is the density of pinning centers, 
C66 = eo/4ao is the vortex-lattice shear modulus; ei = e 2 eo is the 
vortex line tension, e is the anisotropy parameter, eo is related 
to the energy of a vortex line per unit length eo = ($o/4ttA) 2 , 
ao = \f^o/B is the vortex line spacing, $0 is the flux quan- 
tum, A is the penetration depth, £ is the coherence length. The 
characteristic longitudinal length Lo, which determines the size 
of elastic screening of local distortions of the vortex line can be 
evaluated by minimizing of the elastic energy with respect to L 



The experiments of magnetic relaxation are the most exten- 
sively used tools to study the vortex dynamics in a variety of 
superconducting materials [121-124]. The relaxation rate S as 
a function of H can provide relevant information about the pin- 
ning mechanism [121]. It has been demonstrated [121] that the 
relaxation rate drops as field increases from below B 3V up to 
a field close to second magnetization peak field B sp , increasing 
again as H become larger than B 3p . Results of this magnetic re- 
laxation experiments have shown the existence of a crossover in 
the pinning mechanism at B sv , where apparently, this crossover 
occurs without a change in the behavior of the relaxation rate 
with field. The peak in S(H) has been recently interpreted [121] 
as a consequence of the crossover from individual to collective 
regime. At the high-temperature side of the peak in S(H), re- 
laxation is anomalous in the sense that the relaxation rate varies 
non-monotonically with time. It has been shown [121] that the 
anomalous relaxation that occurs on the high-temperature side of 
the peak in S(H) indicates that the vortex structure undergoes a 
change in the relaxation regime. The second magnetization peak 
occurring at B sp is formed by a crossover in the pinning mecha- 
nism, from collective to plastic pinning as the field increases [122]. 
Results of this work show that the temperature dependence of 
the second magnetization peak position, B ap is well explained in 
terms of a plastic motion of the vortex lattice. An analysis of the 
relaxation data shows that the peak field is moving with time to 
lower values. The dynamic picture is also consistent with results 
by Yeshurun et al. [130], showing a clear time dependence of 
the second peak. Both studies show that the peak is virtually 
absent when the loops are measured on a short time scale while 
the peak gradually appears as the time scale is expanded. As 
pointed out in the previous section, the activation energy grows 
with B below the second peak and decreases after the peak. This 
crossover in the field dependence indicates an elastic-to-plastic 
creep crossover around the peak-field similar to that observed in 
other materials [122-130]. 
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The characteristic elastic energy of the vortex thus becomes 
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The collective pinning length can be evaluated by balancing the 
elastic energy cost and the pinning energy gain associated with 
a small displacement of a region of linear size L 
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Equating the resulting pinning force and the Lorentz force on 
segments of the length L c allows one to obtain the critical cur- 
rent density for the case of weakly pinning point defects and 
small applied fields 
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Let us estimate for the single-vortex — small-bundle crossover 
field B sv , where the system changes from single- vortex creep to 
collective creep of small bundles. This crossover field B sv is de- 
termined by the condition [49] 

L c = ea . 
This condition is equivalent to 



B sb = 



Jo 



-B C 2, 



(108) 



where /3 « 5 is the constant parameter, jo = 4_B c /3v / 6/^oA the 
depairing current density, B c = $o/2v27rA£ the thermodynamic 
critical field, B C 2 = pio$o/27r£ 2 the upper critical field. 

As can be seen from the last expression, the temperature de- 
pendence of the crossover magnetic field B sv is determined by the 
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temperature dependence of the intrinsic superconducting param- 
eters, such as the coherence length £(T) and the magnetic pene- 
tration depth A(T), and the disorder parameter 7 of the vortex 
system. According to G. Blatter et al. [49], defects can interact 
with the vortices in two different ways. They may cause a spa- 
tial variation of the transition temperature (<5T c -pinning) and a 
spatial modulation of the mean free path (<5/-pinning). In both 
cases, the influence of disorder is described by a disorder param- 
eter 7, proportional to the defect density and the temperature 
dependence of 7 is different for the two cases. For 5T c -pinning, 
7 = A~ 4 , while for <5Z-pinning, 7 = A -4 . Thus, depending on the 
type of pinning the crossover field can be either an ascending or 
descending function of temperature. Considering the functional 
form of the disorder parameter for 5T c -pinning, 7 = A" 4 . Within 
the frame of the Ginzburg-Landau theory A(T) = (1 — t 4 )~ 1 ^ 2 
and£(T) = [(l-i 2 )/(l-i 2 )]- 1/2 , where A(0) and £(0) is the pen- 
etration depth and the coherence length at T=0 K. Inserting the 
values of the coherence length f(T) and the penetration depth 
A(T) into the last equation for B ap (T) we obtain the following 
expression for the crossover field B sv in the single creep regime 
for the case 8T C pinning 



B sv = B 3V (0) ( 1 + t2 
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and <5Z-pinning, respectively 



B sv = B sv (0) 



1 - t 
1 + t 2 
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Thus, the vortex system is in the single- vortex pinning regime at 
low fields. 

As we have already mentioned above, the analysis of collec- 
tive creep in terms of individually moving vortices applies only 
for a limited regime at B — B sv . However, with increasing ap- 
plied field the relative importance of the interaction between the 
vortices grows and for larger field values, B sv < B we enter to 
the collective or elastic creep, where thermal fluctuations become 
large due to large flux creep [122-130]. It was now strongly ex- 
perimentally established that, the line B ap can be well explained 
by a plastic creep model based on dislocation mediated motion of 
vortices similar to diffusion of dislocations in atomic solids [122- 
130]. The activation energy U p i for the motion of a dislocation 
in the vortex lattice can be estimated as 



U pl = 



1 

4-7T7 



eodo oc B 



-1/2 



(111) 



One notices that U p i decreases with field in contrast to the collec- 
tive creep activation energies U e i, which increases with field. The 
fishtail peak location B sp can be determined from the condition 
U p i — Uei- Using the logarithmic solution of the flux diffusion 
equation we find 



kT In- 
to 



(112) 



Balance between the U p i and £/ e ; energies gives an expression for 
the temperature dependence of B sp . To roughly estimate the 
crossover field and equate the values of the elastic and plastic 
pinning barriers [122], the following relation can be used 



f'o 



5 T c - pinning 
5 1 - pinning 
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Fig. 12. Crossover field as a function of temperature 
corresponding to 51 and 8T C pinning. 

The data in Figure 12 clearly show that only, the 5T c -pinning 
formula can explain the temperature dependence onset field B on , 
accurately. Thus, a crossover from individually pinned regime to 
a collectively pinned regime is to be anticipated at a threshold 
field B sv , above which vortex- vortex interactions become domi- 
nant. At low temperatures with the numerical values A and £ we 
find that B sv is of the order of 1.7 T. This is justified in the lim- 
iting case of small magnetic fields. Consequently, the transverse 
correlation length of the flux lattice L c =l is approximately equal 
to O.lao, i.e. it is much lower than the intervortex spacing oq. 



= kT\nt. 



From the last relation we find 



B sp = $0 



„2„2 

e e 
{kT) 2 



2 

oc e . 



(113) 



(114) 



Inserting the values of the penetration depth A to the equation 
for the B 3p we obtain the temperature dependence of the transi- 
tion field, as it is demonstrated in Fig. 13. Note, that a similar 
temperature dependence B sp (T) was found in some experiments 
[116] for the melting line. This interesting similarity may sup- 
port previous claims [117] that the plastic motion of defects in 
the vortex lattice is a precursor to the melting transition. 
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Fig. 13. Second magnetization peak field as a function of 
temperature. 



To describe the order-disorder phase transition line, we use the 
following Lindemann criterion 



In Fig. 13 we show the temperature dependence of the peak 
field. Thus, our analysis of the vortex-dynamics shows that 
the onset of the second magnetization peak occurring at B on 
is formed by a crossover in the pinning mechanism, from single 
to collective pinning as the field increases. On the other hand, 
the temperature dependence of the second magnetization peak 
B sp can be well explained by the last expression for the plastic 
vortex creep model. 



§5. 3. Order disorder transition 

The nature of the different vortex phases in the mixed state 
of high temperature superconductors is a topic of extensive ex- 
perimental and theoretical research [49, 143-146]. Much of the 
experimental work has been focused on the highly anisotropic 
(BSCCO) crystals, revealing a rich phase diagram [143, 145, 146]. 
In particular, recent local magnetization measurements [143] in 
BSCCO crystals revealed a sharp onset of an anomalous sec- 
ond magnetization peak at a field B on , which was interpreted 
as indicating a transition between two solid phases of the vortex 
structure. An evidence for two distinct solid vortex phases in 
BSCCO was previously obtained in neutron diffraction [145] and 
mSR [146] experiments. Following these observations, a theoret- 
ical model was developed [147-149], describing a mechanism for 
a disorder-induced transition, from a relatively ordered vortex 
lattice, to a highly disordered vortex state. 

An experimental studies have shown [153, 154] that the sam- 
ple exhibits a distinct second peak, i.e., a strong increase in the 
magnitude of the magnetization in an intermediate field range. 
The lowest of these is B on denoting the onset field of the sec- 
ond peak while the peak maximum occurs at B 3p . At low fields 
the elastic interactions govern the structure of the vortex solid, 
forming a quasiordered lattice [150]. Above the onset field B on , 
disorder dominates and vortex interactions with pinning centers 
result in an entangled solid where cells of the vortex lattice are 
twisted and dislocations proliferate [147-149]. To understand the 
origin of the onset peak it is very useful to study the tempera- 
ture dependence of the field B on . By considering the compe- 
tition between the elastic energy of the vortex system and the 
pinning energy [147-149] it has been suggested that the second 
peak can result from a transition of a low-field quasiordered vor- 
tex phase to a disordered vortex solid at higher fields, induced 
by the quenched disorder. According to recent theoretical model 
[147-149], the vortex phase diagram is determined by the in- 
terplay between three energy scales: the vortex elastic energy 
E e i, the energy of thermal fluctuations E t h and the pinning en- 
ergy E p . The competition between the first two determines the 
melting line while the competition between the last two deter- 
mines the irreversibility line. The competition between the elas- 
tic energy and the pinning energy determines the order-disorder 
transition field B on . We now study the temperature dependence 
of the onset field, within the framework of a theoretical model 
[147-149], based on the phenomenological Lindemann criterion 
for the regime of single vortex pinning [52]. The order-disorder 
occurs when the disorder-induced wandering of the flux lines is 
comparable to the lattice constant do and the vortex system loses 
its translational order and transforms into an disordered phase, 
in which vortices better adapt to the local pinning potential. 



<u 2 (L ,0)> = c\al, (115) 

where (u) describes the mean relative displacement of neighbor- 
ing flux lines caused by the disorder. This Lindemann crite- 
rion leads to the usual condition for the order-disorder transition 
at low temperatures, when the thermal fluctuation is negligible. 
The Lindemann number cl is introduced here as a phenomeno- 
logical parameter that is supposed to depend only weakly on the 
specific lattice parameters of the solid phase, in particular it is 
assumed to be independent of the magnetic field. Strictly speak- 
ing, the value of the constant cl may depend on whether the 
order-disorder transition occurs in the single vortex pinning re- 
gion or in the region of bundle pinning. However, to understand 
the essence of the matter, we shall use the simplest approxima- 
tion: cl will be considered as the same constant for the various 
regimes of pinning. We, therefore, use this value in our calcula- 
tions. 

The physics of a single vortex line in point disorder exhibits 
two different scaling regimes depending on the typical size of the 
disorder-induced mean-square displacement 



(ML) 



■u(0)] 2 )=u 2 



1 2 C 



(116) 



of a vortex segment of length L, where < > denotes the 

full statistical average to be taken over dynamical variables first 
and then over the quenched variables describing the disorder, £ 
is so-called static critical exponent. We have to carefully dis- 
tinguish several pinning regimes depending on the size of the 
pinning length L c in comparison to the single- vortex length Lq. 
For short vortex lengths Lq < L c displacements are small 

(u 2 (L c )) ~ e, 

perturbation theory is valid, and we can work with Larkin's ran- 
dom forces [141]. Fluctuations around the ground state of the 
line are Gaussian in this random force regime and we find a 
roughness exponent z = 3/2, i.e., 



([u(L)-u(0)] 2 )=e 



This regime is valid up to the collective pinning or Larkin length 
L c which is defined by the condition 

Longer segments explore many almost degenerate minima of 
the pinning energy landscape such that fluctuations are non- 
Gaussian. In this so-called random manifold regime the rough- 
ness exponent is known exactly z — 5/2 [155]. The currently 
most reliable estimate for general n is z = 5/8 for the vortex line 
in the bulk superconductor. Therefore we find 



<[u(L) - u(0)} 2 ) = ui 



L_ 



-,5/4 



On scales exceeding the Larkin length L c < Lq, the vortices 
start to explore many minima of the disorder potential [155]. 
This regime is referred to as the random manifold regime [155, 
151]. In the case L c < Lq the pinning energy becomes 
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where 



E p = U, 



Ud P — 



d p 



1 1/5 



(117) 
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By comparing the pinning E p and elastic energies E e i we now 
derive the transition field 



B = 



B on = Bo 



T 



U, 



dp 



(118) 



To = 



ceeoC 



^2 > 2 

Each vortex remains individually pinned in the presence of ther- 
mal fluctuations until the thermal energy T is greater than the 
typical depinning energy Ud P of each single vortex. After further 
analytical calculation we obtain the following formula to describe 
the temperature dependence of the transition field 



Bo 



(119) 



As we have mentioned above, depending on the type of pinning 
the crossover field can be either an ascending or descending func- 
tion of temperature. Substituting the temperature dependencies 
of the quantities £(T), A(T) into the relation (119) and consid- 
ering the functional form of the disorder parameter, we arrive at 
the following expressions for <5T c -pinning 



B on oc (1 — t ) 



4x3/2 



and for <5/-pinning 



respectively. 
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(120) 
(121) 



It is notice able, that onset field B on is inversely proportional 
to both the anisotropy 7 and the disorder parameter 7. Thus, 
both the anisotropy and disorder may sufficiently affect to the 
shape of temperature dependencies of both the onset and second 
magnetization peak field. As has been pointed out Yeshurun 
et al. [130] that the point disorder induced by electron irra- 
diation can modify the penetration depth A, anisotropy 7 and 
critical temperature T c , and thus cause a significant shift of the 
transition line. Thus, the line B n(T) can be expected to shift 
lower fields with the introduction of additional disorder in the 
crystal structure. Calculations [147-149] of B on , based on a Lin- 
demann criterion, confirm this expectation. Other experimental 
studies Khaikovich et al. [143] of the effect of electron irradia- 
tion on the order-disorder transition line have found a system- 
atic decrease of _B on (T) with increasing irradiation dose. The 
decrease in B on is consistent with the enhanced vortex pinning 
after electron irradiation. Therefore, the quasiordered phase is 
stabilized by introduction of the point disorder and the strong 
pinning region is expanded to the lower fields. These results of 
electron irradiation effect provide further evidence of the field- 
driven disordering transition scenario [147-149] as a possible ori- 
gin of B n(T)- Specifically, the decrease in the superfluid den- 
sity n a with increasing doping at high temperatures increases 
the magnetic penetration length A, hence decreases B on oc A~ 4 . 
It can be concluded that introduction of even a small amount 
of disorder into very clean systems changes the phase diagram 
drastically by lowering the order-disorder transition line B on (T) 
and also changing it's temperature dependence. When the disor- 
der reaches a certain threshold, however, introducing additional 
disorder does not continue this tendency. With other words, in 
highly disordered systems, the vortex matter phase diagram is 
relatively robust with respect to variations in the exact degree 
of disorder [156-158]. 

Conclusion 
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In the present paper we studied the flux jump instabilities of the 
critical state in conventional and high-T c superconductors. We 
have determined the flux jump field and critical state stability 
criterion within framework of Bean's model. To determine flux 
jumps threshold an analytical simulations of coupled equations 
for the magnetic, electric field inductions and temperature has 
been performed. The field of the first flux jumps is calculated an- 
alytically using the dynamic and the adiabatic approximations. 
The qualitative analysis of the magnetic flux jumps instabilities 
for superconductors is provided. 
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Fig. 14. Onset field as a function of temperature. 

In Fig. 14 we show the temperature dependence of the onset 
field. The data in Figure 14 clearly shows that B n(T) for 5T C 
pinning is a decreasing functions of temperature and therefore 
this formula can explain the onset of the second magnetization 
peak, properly. 
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